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Introduction: 


This  project  aims  to  study  interactions  between  genetic  and  environmental  factors  in  a 
viable  system  -  human  fibroblasts.  Fibroblasts  with  a  scleroderma  (SSc)  susceptible  genetic 
background  may  be  more  vulnerable  to  environmental  triggers.  Studies  of  fibroblasts  with  and 
without  SSc  susceptible  backgrounds  in  response  to  potential  environmental  triggers  will  provide  a 
great  opportunity  to  understand  etiopathogenesis  of  SSc. 

Body: 

According  to  our  SOW,  there  are  four  major  tasks  that  will  be  performed  for  this  project.  1.  We 
will  culture  and  genotype  human  fibroblasts  obtained  from  skin  biopsies.  2.  We  will  stimulate 
cultured  fibroblasts  with  silica,  cellulose,  titanium  oxide  and  carbon  particles  individually,  and  then 
examine  potential  alterations  of  fibroblasts  in  terms  of  proliferation,  ROS  production,  collagen 
synthesis,  cytokine  expression  and  release.  3.  We  will  apply  comprehensive  statistical  analysis  to 
identify  the  association  of  genetic  susceptibility  genes  and  these  environmental  factors  in  terms  of 
the  degree  of  their  profibrotic  effects.  4.  We  will  also  study  specific  biological  pathways  triggered  by 
risk  elements  in  cultured  fibroblasts. 

During  this  second  year  of  the  project,  we  have  been  continuously  obtaining  skin  biopsies  from 
SSc  patients  and  normal  controls.  Currently,  we  have  a  total  of  166  skin  fibroblast  strains  including 
65  from  SSc  patients  and  101  from  normal  controls.  These  primary  fibroblast  cultures  were 
established  from  each  skin  biopsy,  and  maintained  in  healthy  condition  for  proposed  functional 
studies. 

We  have  performed  stimulation  assays  on  82  primary  fibroblast  strains  for  part  of  our  proposed 
studies.  Each  of  these  82  fibroblast  strain  was  stimulated  with  silica  particles  for  evaluation  of  their 
responses  toward  potentially  fibrotic  changes  over  5  different  time-points.  Meanwhile,  we  have 
genotyped  79  of  these  fibroblast  strains  with  both  the  HLA  subtyping  and  the  genome-wide  single 
nucleotide  polymorphism  (SNP)  profiling  (lllumina  Human  317  K  SNP  panel).  It  should  be  noted 
that  the  studies  of  performing  genome-wide  SNP  profiling  were  funded  by  our  other  research 
projects.  This  addition  of  genotyping  appeared  to  be  very  important  in  understanding  complexity  of 
human  genetics,  especially  in  association  with  environmental  factors.  We  examined  gene 
expression  changes  of  COL1A2,  COL3A1,  CTGF,  SPARC  and  TIMP3  of  these  fibroblasts  in  five 
different  time  points  (0,  1,  2,  3,  and  5  days)  in  response  to  silica  stimulation.  The  COL1A2, 
COL3A1,  CTGF,  SPARC  and  TIMP3  are  three  major  components  of  the  extracellular  matrix 
(ECM),  which  are  actively  involved  in  fibrotic  changes  of  the  tissues.  Using  longitudinal  linear 
models  in  the  analysis  of  association  between  specific  genotypes  and  dynamic  changes  of 
gene  expression  of  the  fibroblasts  in  responding  to  silica  stimulation,  we  identified  that 
specific  SNPs  were  associated  with  either  single  gene  expression,  or  paired  gene 
expression  (Table  1).  Notably,  the  paired  gene  expressions  of  CTGF/SPARC  and 
COL1A2/COL3A1  were  associated  with  multiple  SNPs,  which  suggest  a  strong  biological 
correlation  between  these  genes.  HLA-DQB1*0301  is  a  positive  marker  for  SSc 
susceptibility.  Interestingly,  HLA-DQB1*0301  positive  fibroblast  strains  showed  a 
significant  up-regulation  (>2-fold  increase  of  transcript)  of  COL1A2,  COL3A1,  CTGF  and 
SPARC  in  response  to  silica  stimulation  after  24  hours,  which  was  significantly  different 
from  the  response  of  HLA-DQB1*0301  negative  fibroblast  strains  (P  <  0.05  for  all  4 
genes). 

While  HLA-DQB  1*0301  is  a  previously  identified  SSc  genetic  marker,  we  are  recently 
exploring  novel  SSc  susceptibility  markers  using  a  genome-wide  SNP  association  (GWAS) 
approach.  Collaborating  with  Dr.  Eun  Bong  Lee  (Seoul  National  University  College  of 
Medicine,  Seoul,  Korea),  we  demonstrated  that  specific  SNPs  of  the  HLA-DPB1  and  -DPB2  were 
strongly  associated  with  SSc  susceptibility  in  both  Koreans  and  US  Caucasians  (see  attached 
manuscript  1).  Incorporating  this  information  into  our  data  of  gene  expression  changes  of  the 
fibroblasts,  we  revealed  that  specific  SNPs  of  HLA-DPB1  also  was  associated  with  gene 
expression  changes  of  the  COL1A2  of  the  fibroblasts  in  response  to  silica  stimulation  (p  = 
0.00042).  Moreover,  several  other  potential  SSc  susceptibility  regions  identified  from  the  Dr.  Lee’s 
GWAS  also  showed  an  overlap  with  the  genetic  regions  associated  with  gene  expression  changes 


of  the  ECM  components  of  the  fibroblasts  in  responds  to  silica  stimulation.  The  figure  1 
demonstrates  an  example  of  the  association  between  genetic  regions  and  expression  changes  of 
the  COL1A2  of  the  fibroblasts  in  response  to  silica  stimulation.  The  figure  2  represents  a  GWAS 
results  from  Dr.  Lee’s  data  that  shows  log  P-values  of  the  genome-wide  SNPs  in  association  with 
SSc.  Comparing  both  figures,  we  can  find  the  overlapped  chromosomal  regions  such  as  on 
chromosome  1,  3,  6  and  10.  In  addition  to  Dr.  Lee’s  GWAS  studies  of  SSc,  Dr.  Maureen  Mayes,  a 
professor  in  our  division  of  Rheumatology  in  UT  Health  Science  Center,  just  completed  the  GWAS 
studies  of  SSc  with  1678  SSc  Caucasian  SSc  patients  and  5530  Caucasian  controls  (Dr.  Mayes’ 
personal  communication).  We  had  a  chance  to  pre-examine  this  unpublished  GWAS  data. 
Remarkably,  multiple  SNPs  that  showed  an  association  with  SSc  susceptibility  in  the  GWAS 
studies  were  overlapped  with  the  SNPs  that  also  were  associated  with  the  gene  expression  of  the 
ECM  of  the  fibroblasts  in  response  to  silica  stimulation  (Table  2).  Further  studies  of  the  genes 
corresponding  to  these  SNPs  may  provide  important  knowledge  to  understand  pathogenesis  of 
SSc.  For  instance,  the  ABCA7  (ATP-binding  cassette,  subfamily  A,  member  7)  encodes  a 
transmembrane  protein  involved  in  energy-dependent  transport  of  a  wide  spectrum  of  substrates. 
Dysregulation  of  ATP-dependent  membrane  transportation  may  lead  to  various  pathological  events 
associated  with  vascular  and  immune  functions.  The  APBA1  (amyloid  beta  A4  precursor  protein¬ 
binding,  family  A,  member  1)  is  a  candidate  gene  for  Alzheimer  disease.  However,  serum  amyloid 
protein  was  reported  to  be  lower  in  SSc  patients  ( Pilling  D.,  et  al.  J  Immunol.  2003;  1 7 1:5537- 
46).  According  to  Philing’s  studies,  low  levels  of  SAP  may  augment  pathological  processes 
leading  to  fibrosis.  The  EGKD  (diacylglycerol  kinase  delta)  is  a  gene  encoding  protein 
regulates  protein  kinase  C  and  epidermal  growth  factor  receptor  signaling  (Proc.  Nat.  Acad. 
Sci.  103:  15485-15490,  2006/  The  RBMS3  (RNA-binding  motif  protein,  single  strand¬ 
interacting  3)  is  a  gene  closely  related  to  c-myc  involving  in  regulation  of  DNA  replication, 
gene  transcription,  apoptosis  and  cell  cycle. 

Overall,  above  data  strongly  support  our  original  proposal  that  genetic  elements  within  SSc 
fibroblasts  might  contribute  to  susceptibility  to  fibrotic  process.  Integrative  studies  of  genetic  and 
environmental  factors  with  human  fibroblasts  may  facilitate  the  discovery  of  potential  pathogenesis 
of  SSc. 

Our  efforts  in  studying  specific  biological  pathways  triggered  by  risk  elements  in  cultured 
fibroblasts  have  been  continued  through  last  year.  Previously,  we  demonstrated  that  silica 
stimulation  on  fibroblasts  induced  TGF-(3  signaling  (reported  in  original  proposal),  which  may 
contribute  to  fibrosis.  Recently,  we  identified  that  SPARC  inhibition  attenuated  profibrotic  effects  of 
TGF-p  on  human  fibroblasts.  Last  year,  we  applied  SPARC  siRNA  in  bleomycin-induced  mouse 
fibrotic  model,  and  demonstrated  an  anti-fibrotic  effect  of  SPARC  inhibition  in  vivo  (attached 
manuscript  2).  Although,  the  in  vivo  studies  of  anti-fibrosis  were  funded  by  other  projects,  the 
results  represent  a  successful  example  of  translational  studies  of  medical  diseases. 


Table  1.  Associations  between  SNPs  and  gene  expression  of  C0L1A2,  C0L3A1.  CTGF, 
SPARC  and  TIMP3  of  human  fibroblasts  in  response  to  silica  stimulation 


P-value  fortesting  genetic  effect  to  gene  expression 


Associated  genes 

rs  number 

Gene  Name 

COL1A2 

COL3A1 

CTGF 

SPARC 

TIMP3 

COL1A2/COL3A1 

rsl0484710 

COL21A1 

6  XlO'12 

1.7  XlO'5 

>0.01 

>0.01 

>0.01 

COL1A2/COL3A1 

rs5972761 

DMD 

7.4  xlO'12 

8.8  XlO-4 

>0.01 

>0.01 

>0.01 

COL1A2/COL3A1 

rs5972763 

DMD 

7.4  XlO'12 

8.8  XlO-4 

>0.01 

>0.01 

>0.01 

COL1A2/COL3A1 

rs3734061 

FAT2 

5.4  XlO'12 

7.9  XlO-4 

>0.01 

>0.01 

>0.01 

COL1A2/COL3A1 

rs998076 

FAT2 

5.4  xlO-5 

3  XlO-4 

>0.01 

>0.01 

>0.01 

COL1A2/COL3A1 

rs2288777 

FAT2 

3.7  xlO'4 

1.3  xlO'5 

>0.01 

>0.01 

>0.01 

COL1A2/COL3A1 

rsl410871 

GBP4 

7.1  XlO'8 

6.5  XlO"4 

>0.01 

>0.01 

>0.01 

COL1A2/COL3A1 

rs2196426 

hCG_2015138 

7.2  XlO14 

1  XlO'7 

>0.01 

>0.01 

>0.01 

COL1A2/COL3A1 

rs4803831 

OPA3 

4.2  XlO'8 

1.3  XlO'5 

>0.01 

>0.01 

>0.01 

COL1A2/COL3A1 

rs3809251 

PARP11 

7.9  XlO'8 

1.3  XlO'7 

>0.01 

>0.01 

>0.01 

COL1A2/COL3A1 

rs2061783 

XRCC4 

2.9  XlO'14 

6.3  XlO'8 

>0.01 

>0.01 

>0.01 

COL1A2/COL3A1 

rsl0947541 

DEF6 

2.9  XlO'14 

6.3  XlO'8 

>0.01 

>0.01 

>0.01 

COL1A2/COL3A1/SPARC 

rs2289722 

C4orf20 

2.1  XlO"5 

3.6  XlO'12 

>0.01 

7.5  XlO"5 

>0.01 

CTGF/SPARC 

rs6810854 

ACOX3 

>0.01 

>0.01 

3.4  XlO'15 

8.8  XlO-6 

>0.01 

CTGF/SPARC 

rs9375513 

C6orfl74 

>0.01 

>0.01 

1.4  xlO'8 

1.4  XlO'5 

>0.01 

CTGF/SPARC 

rs2236026 

KIAA0408 

>0.01 

>0.01 

1.4  xlO'8 

1.4  xlO'2 

>0.01 

CTGF/SPARC 

rs3922703 

PLCXD2 

>0.01 

>0.01 

1.5  XlO'8 

9.4  XlO'5 

>0.01 

CTGF/COL3A1 

rs3788926 

ZRSR2 

>0.01 

1.1  XlO'8 

2.8  XlO'5 

>0.01 

>0.01 

CTGF/TIMP3 

rsll702035 

C21orf51 

>0.01 

>0.01 

1.4  XlO"7 

>0.01 

8.8  XlO'5 

CTGF/SPARC/TIMP3 

rsl0507399 

HSPH1 

>0.01 

>0.01 

7.8  XlO"9 

2.9  XlO-6 

3.9  XlO'8 

SPARC/CTGF/COL1A2/COL3A1 

rsl3082485 

PIK3CA 

2.5  XlO'7 

1.8  XlO'2 

2.5  XlO'2 

1.9  XlO'7 

>0.01 

COL1A2 

rs437444 

CELSR2 

1.1  XlO'7 

>0.01 

>0.01 

>0.01 

>0.01 

COL1A2 

rsl0478113 

MCC 

1  XlO'7 

>0.01 

>0.01 

>0.01 

>0.01 

COL1A2 

rs3793726 

PRKCQ 

1.5  XlO'7 

>0.01 

>0.01 

>0.01 

>0.01 

COL1A2 

rs7918923 

PRKCQ 

1.6  XlO'7 

>0.01 

>0.01 

>0.01 

>0.01 

COL1A2 

rs6445245 

PTPRG 

1.5  XlO'7 

>0.01 

>0.01 

>0.01 

>0.01 

COL1A2 

rs2864780 

RAB31 

7.8  XlO'8 

>0.01 

>0.01 

>0.01 

>0.01 

COL1A2 

rs3731895 

RESP18 

6  XlO'8 

>0.01 

>0.01 

>0.01 

>0.01 

COL3A1 

rsllS 64394 

CDH2 

>0.01 

1.5  XlO"7 

>0.01 

>0.01 

>0.01 

COL3A1 

rsll754507 

DSE 

>0.01 

1.6  XlO'11 

>0.01 

>0.01 

>0.01 

COL3A1 

rs 70529 34 

SPANXN3 

>0.01 

3.4  XlO'8 

>0.01 

>0.01 

>0.01 

COL3A1 

rsll817964 

ZNF365 

>0.01 

1.6  XlO'7 

>0.01 

>0.01 

>0.01 

TIMP3 

rsl0499156 

LAMA2 

>0.01 

>0.01 

>0.01 

>0.01 

5.4  xlO'8 

TIMP3 

rsl3154825 

SUT3 

>0.01 

>0.01 

>0.01 

>0.01 

8.1  XlO'8 

TIMP3 

rs!476707 

C19orf29 

>0.01 

>0.01 

>0.01 

>0.01 

6.1  XlO'8 

Table  2.  The  SNPs  associated  with  both  SSc  susceptibility  and  gene  expression  changes 
of  the  fibroblasts  in  response  to  silica  stimulation. 


SNP  ID 

Gene  symbo 

Gene  ID 

P-value  forSNP- 
SSc  association 

P-value  fortesting 
genetic  effect  on 
gene  expression 

Associated 

gene 

expression  of 

rs3752228 

ABCA7 

NM_019112.2 

0.0004622 

0.0046 

COL1A2 

rs1831555 

APBA1 

NM_001 163.2 

0.0001345 

0.0077 

CTGF 

rs1977552 

APBA1 

NM_001 163.2 

0.0001478 

0.0068 

CTGF 

rs2781530 

APBA1 

NM_001 163.2 

0.0002292 

0.0077 

CTGF 

rs673576 

C21orf88 

NM_1 53754.1 

0.0002405 

0.0057 

TIMP3 

rs7584554 

DGKD 

NM_1 52879.2 

0.0002883 

0.0017 

CTGF 

rs4981200 

NPAS3 

NM_022123.1 

0.0006308 

0.0067 

CTGF 

rs4981200 

NPAS3 

NM_022123.1 

0.0006308 

0.0049 

SPARC 

rs4981200 

NPAS3 

NM_022123.1 

0.0006308 

0.0068 

COL3A2 

rs3773046 

RBMS3 

NM  001003792 

0.0001188 

0.0072 

COL1A2 

Figure  1 .  Identification  of  genetic  loci  associated  with  high  sensitivity  of  fibroblasts  over¬ 
expressing  the  COL1A2  in  response  to  silica  stimulation.  (X  axis  indicates  SNP  positions 
on  different  chromosomes,  Y  axis  indicates  log  p  values  of  association  between  over¬ 
expression  of  the  COL1A2  and  SNP  genotypes).  The  gene  expression  of  the  COL1A2  was 
examined  with  real-time  RT-PCR.  The  genome-wide  SNP  typing  on  human  fibroblasts  was 
performed  with  lllumina’s  human  317  K  SNP  panel. 
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Figure  2.  Identification  of  the  major  loci  associated  with  systemic  sclerosis  in  genome  wide 
scan.  A  total  of  440,734  SNPs  were  evaluated  in  133  patients  with  SSc  and  557  healthy 
controls.  Distribution  of  -Logio  P  values  are  plotted  against  chromosomes.  The  genome¬ 
wide  SNP  typing  was  performed  with  Affymetrix  Whole  Genome  Human  SNP  Array  5.0. 


Key  Research  Accomplishments 

•  Obtained  a  total  of  166  human  fibroblast  strains  (65  SSc  patients  and  101  normal  controls). 

•  Completed  silica  stimulation  in  82  fibroblast  strains  and  obtained  RNA  and  protein  extracts 
from  each  of  the  experiments. 

•  Completed  genotyping  for  79  fibroblast  strains  (34  SSc  and  45  controls)  with  both  HLA 
subtyping  and  genome-wide  SNP  profiling. 


•  Identified  that  the  paired  gene  expressions  of  CTGF/SPARC  and  COL1A2/COL3A1 
were  associated  with  multiple  SNPs,  which  suggest  a  strong  biological  correlation 
between  these  genes. 

•  Identified  some  novel  genes  and  SNPs  associated  with  both  SSc  susceptibility  and  human 
fibroblast  response  to  silica  stimulation  in  terms  of  their  gene  expression  changes  of  the 
ECM  components. 

•  Demonstrated  that  specific  inhibition  of  SPARC,  a  major  ECM  component,  attenuated 
mouse  fibrosis  induced  by  bleomycin  in  vivo. 

Reportable  Outcomes 

During  this  second  year  period,  we  reported  that  SPARC  inhibition  attenuated  profibrotic 
effects  of  TGF-p  in  vitro  and  bleomycin  in  vivo  (oral  presentation)  in  annual  meeting  of  American 
College  of  Rheumatology  2008.  We  have  one  full  article  accepted  by  the  European  Journal  of 
Human  Genetics  (2009).  We  have  two  manuscripts  submitted  for  publication  (one  received 
comments  for  revision)  and  two  manuscripts  in  preparation. 

Importantly,  we  established  a  total  of  166  human  fibroblast  strains  (65  SSc  and  101 
controls),  which  is  close  to  our  minimum  target  of  200  (100  SSc  and  100  controls)  that  can  be  used 
for  many  studies  toward  human  health. 
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Space  Model  in  Studies  of  human  Complex  Disease.  PLoS  ONE,  2008  Feb, 
3(2):e1693. 

2.  Momiao  Xiong,  Gang  Peng  ,  Li  Luo  ,  Yun  Zhu  ,  Pengfei  Hu  ,  Shengjun  Hong  ,  Jinying 
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2.  JC  Wang,  Xiong  H,  Tan  FK,  Arnett  FC,  Zhou  XD.  Studies  of  genetic  and 
environmental  factors  in  human  fibroblasts  showed  novel  gene-gene  interactions. 


Conclusion 

During  this  second  year  of  the  project,  we  established  multiple  primary  cell  strains  from 
normal  controls  and  SSc  patients.  We  performed  stimulation  assays  with  silica  in  82 
primary  fibroblast  strains.  Our  results  showed  that  silica  activate  fibroblasts  toward  fibrotic 
changes.  However,  different  fibroblast  strains  obtained  from  different  individuals  showed 
different  responses  in  terms  of  the  gene  expression  of  the  ECM  components.  Using 
longitudinal  linear  models  in  analysis  of  association  between  specific  genotypes  and 
dynamic  changes  of  gene  expression  of  the  fibroblasts  in  responding  to  silica  stimulation, 
we  identified  that  specific  SNPs  were  associated  with  either  single  gene  expression,  or 
paired  gene  expression  such  as  CTGF/SPARC  and  COL1A2/COL3A1,  which  suggest  a 
strong  biological  correlation  between  these  genes.  Moreover,  some  SNPs  and/or  their 
corresponding  genes  were  found  to  be  associated  with  both  SSc  susceptibility  and  the 
fibrotic  changes  of  human  fibroblast  in  response  to  silica  stimulation.  These  SSc 
susceptibility  genes  include  not  only  previously  identified  ones,  but  also  some  novel  ones, 
such  as  HLA-DPB1  and  APBA1  .These  observations  supported  our  original  proposal  that 
genetic  elements  within  SSc  fibroblasts  might  contribute  to  susceptibility  to  fibrotic  process. 
Integrative  studies  of  genetic  and  environmental  factors  with  human  fibroblasts  may  facilitate  the 
discovery  of  potential  pathogenesis  of  SSc.  In  this  year,  we  will  continuously  obtain  more  human 
fibroblasts,  especially  SSc  fibroblasts.  We  also  will  continuously  perform  stimulation  assays  of 
newly  obtained  human  fibroblasts  for  a  better  chance  to  identify  genetic  components  inside  SSc 
patients  contributing  to  susceptibility  to  environmental  stimuli.  Meanwhile,  we  will  try  to  explore 
which  specific  bio-pathways  associated  with  which  specific  genetic  factors  involved  in  silica 
induced  fibrotic  changes.  Therefore,  our  studies  are  fulfilled  with  original  proposal  in  the  grant. 


Appendices  (two  manuscripts  submitted  for  publication) 
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Abstract 


We  performed  a  genome-wide  association  study  using  the  Affymetrix  Human 
SNP  Array  5.0  in  137  patients  with  systemic  sclerosis  and  564  controls  from 
Korea.  The  SNPs  (rs31 28930,  rs7763822,  rs7764491,  rs31 17230  and 

rs31 28965)  of  HLA-DPB1  and  -DPB2  on  chromosome  6  formed  a  distinctive 
peak  with  log  p-values  (p  <  1  x  10'14)  for  association  with  SSc  susceptibility. 
Subtyping  analysis  of  HLA-DPB1  showed  that  DPB1*1301  (p  =  7.61x10'8)  and 
DPB1*0901  (p  =  2.56x1  O'5)  are  the  most  susceptible  subtypes  for  SSc  in  Koreans. 
We  then  replicated  the  results  in  1,106  SSc  patients  and  2,747  controls  from  a 
US  Caucasian  population.  Remarkably,  two  pairs  of  SNPs,  rs7763822/rs7764491 
and  rs31 17230/rs31 28965,  showed  strong  association  with  Caucasian  SSc 
patients  who  have  either  circulating  anti-DNA  topoisomerase  I  or  anti-centromere 
autoantibodies,  respectively.  Therefore,  our  GWAS  and  confirmatory  studies 
revealed  that  the  region  of  HLA-DPBf  and  -DPB2  contains  the  most  susceptible 
loci  to  SSc,  especially  those  patients  with  anti-topoisomerase  I  or  anti¬ 
centromere  autoantibodies. 
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Introduction 


Systemic  sclerosis  (SSc)  is  a  rare  and  complex  connective  tissue  disease 
of  unknown  etiology  characterized  by  fibrosis  and  vasculopathy  of  skin  and 
internal  organs,  as  well  as  several,  mutually  exclusive,  disease-  specific 
circulating  autoantibodies1.  SSc  can  be  clinically  sub-classified  based  on  patterns 
of  skin  fibrosis  into  limited  and  diffuse  forms.  In  addition,  the  majority  of  SSc 
patients  (90%)  have  circulating  anti-nuclear  autoantibodies  (ANA)2.  The  three 
most  common  autoantibodies  (auto-Abs)  are  anti-DNA  topoisomerase  I  (topo  I), 
anti-RNA  polymerase  III,  and  anti-centromere  antibodies,  in  which  the  first  two 
auto-Abs  tend  to  be  associated  with  diffuse  SSc2: 3,  the  last  one  being  strongly 
correlated  with  limited  SSc,  although  these  associations  are  not  complete2’4. 

Genetic  predisposition  is  widely  believed  to  contribute  to  SSc.  Several 
genetic  regions  have  been  reported  to  be  associated  with  SSc  susceptibility, 
including  the  major  histocompatibility  complex  (MHC)  or  HLA-class  II5,  protein 
tyrosine  phosphatase  non-receptor  22  (PTPN22)6,  allograft  inflammatory  factor  1 
(AIF1)7,  tumor  necrosis  factor  (TNF)8,  cytotoxic  lymphocyte  antigen-4  (CTLA-4)9, 
transforming  growth  factor  (TGF-|3)10  and  connective  tissue  growth  factor 
(CTGF)11.  However,  the  low  prevalence  of  SSc  (approximately  0.00070- 
0.049%)12;  13  and  clinical/serological  heterogeneity  make  genetic  studies  of  SSc 
difficult  with  some  differing  results  reported  for  the  same  genes  in  different  ethnic 
groups.  To  address  such  difficulties,  we  conducted  a  two-step  genetic 
association  study  in  four  independent  populations  to  identify  the  susceptibility 
markers  for  SSc. 
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Material  and  Methods 


Study  Subjects 

We  examined  4  different  ethnic  populations  (Koreans,  Caucasian  Americans, 
African  Americans  and  Hispanics).  Korean  study  population  was  composed  of 
151  SSc  patients  diagnosed  according  to  the  ACR  preliminary  criteria  for  SSc14. 
All  the  patients  were  enrolled  from  Seoul  National  University  Hospital  between 
January  1998  and  2007.  Genomic  DNA  was  extracted  from  whole  blood  using 
standard  methods.  A  total  of  137  cases  which  passed  the  DNA  quality  check 
were  entered  into  the  GWAS  using  Affymetrix  Genome-Wide  Human  SNP  Array 
5.0.  A  total  of  133  cases  which  showed  >95%  of  call  rates,  were  finally  entered 
into  the  case-control  analysis.  The  mean  age  at  diagnosis  was  42  years  ranging 
from  4  to  67  years.  Mean  duration  of  the  disease  was  10  years  and  the  mean 
time  from  diagnosis  to  blood  sampling  was  5  years. 

The  600  healthy  controls  were  randomly  selected  from  10,000  healthy  Koreans 
belonging  to  Korean  Association  Resource  (KARE)  Project,  based  on  the 
frequency-matching  on  sex  with  the  cases.  The  mean  age  of  the  controls  was 
52.5  years.  The  same  platform  (Affymetrix  Genome-Wide  Human  SNP  Array  5.0) 
was  used  for  the  whole-genome  scan  of  the  controls.  After  excluding  cases  with 
low  call  rate  less  than  95%,  mismatched  sex  and  potential  relatives,  a  total  of  557 
controls  were  finally  entered  into  the  case-control  analysis.  The  institutional 
review  board  of  Seoul  National  University  Hospital  approved  the  study  and  all 
patients  and  controls  provided  written  consents. 
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SSc  patients  of  Caucasians,  African  Americans  and  Hispanics  who  met  the  ACR 
criteria  for  SSc  and  corresponding  controls  (sex  and  ethnically  matched)  were 
enrolled  in  the  Division  of  Rheumatology,  University  of  Texas  Health  Science 
Center  at  Houston  (UTHSC-Houston)6’ 1 .  In  addition,  2,300  Caucasian  controls 
were  selected  from  NIH  data  base  of  Genotype  and  Phenotype  (dbGaP)  found  at 
http://www.ncbi.nlm.nih.gov/gap.  The  study  was  approved  by  the  institutional 
review  boards  of  the  University  of  Texas  Health  Science  Center  at  Houston.  All 
the  individual  patients  and  controls  provided  written  consent. 

GWAS  analysis 

Among  the  500,568  SNPs  present  in  Affymetrix  Whole  Genome  Human  SNP 
Array  5.0,  440,734  SNPs  were  accessible  after  excluding  hidden  SNPs.  Q-Q  plot 
was  obtained  in  the  condition  of  p-value  >  0.0001  for  Hardy-Weinberg  equilibrium 
and  call-rate  >  0.95  (Figure  1).  Finally  cluster  quality  analysis  was  performed, 
which  led  to  349,209  SNPs  which  can  be  analyzed.  The  most  significant  SNPs 
were  determined  to  be  rs31 28930,  re7763822,  rs7764491,  rs31 17230  and 
rs31 28965  and  they  went  into  the  fine  mapping  process. 

Fine  mapping 

For  Koreans,  we  performed  a  fine  mapping  focusing  on  HLA-DPB/  and  -DPB2 
regions  with  137  SSc  cases  and  548  healthy  controls  in  whom  DNA  was 
available.  For  HLA-DPB1,  highly  variable  exon  2  of  the  gene  was  DNA- 
sequenced  to  determine  subtype  of  HLA-DPB1.  For  the  other  region  including 
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HLA-DPB2,  total  22  tag  SNPs  were  selected  with  r2  threshold  of  0.8  and  minor 
allele  frequency  over  5%  in  HapMap  Japanese  panel  data  (release  22)  using  the 
Haploview  version  4.1 15;  16.  In  tag  SNPs,  17  SNPs  which  are  included  in 
Affymetrix  SNP  chip  were  forced  to  be  included. 

For  replication  studies,  TaqMan  Assays  were  used  for  SNPs  rs31 28930, 
re7763822,  rs7764491,  rs31 17230  and  rs31 28965  genotyping  with  an  ABI 
7900HT  Fast  Real-Time  PCR  System  in  Caucasian,  African  American  and 
Hispanic  populations.  Genotyping  results  of  all  five  SNPs  passed  quality  tests  for 
Hardy-Weinberg  (p  >  0.001)  and  calling  rate  (>95%).  Two  groups  of  Caucasian 
controls  from  the  NIH  data  base  and  the  UTHSC-Houston  showed  concordant 
association  with  Caucasian  SSc  patients. 

Statistical  analysis 

The  association  of  specific  SNPs  with  the  disease  or  a  subset  of  the  disease  was 
analyzed  by  the  comparison  of  minor  allele  frequencies  of  the  cases  and 
controls,  with  significance  determined  by  p-values  of  chi-square  tests,  Cochran- 
Armitage  Trend  test,  and  Jonckheere-Terpstra  tests.  The  odds  ratio  of  cases’ 
having  a  selected  SNP  compared  with  the  controls’  and  its  relevant  95% 
confidence  intervals  were  also  determined.  For  GWAS,  The  threshold  for 
declaring  significance  after  Bonferroni  correction  for  adjusting  multiple  tests  was 
p<1.43x10'7.  It  was  p<2.2x10'3  for  the  fine  mapping  study  (n=22)  and  p<0.01  for 
the  replication  study  (n=5).  All  the  association  tests  were  based  on  the 
comparison  of  alleles.  PUNK17  and  SAS  9.1.3  (SAS  Institute  Inc.,  Cary,  NC, 
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USA)  were  used  for  the  statistical  analysis.  Linkage  disequilibrium  analysis  for 
HLA-DPB1  and  -DPB2  regions  was  performed  with  Haploview,  version  4.1 16. 

Results 

We  first  examined  a  Korean  population  who  are  relatively  homogenous  in 
which  62.2%  of  patients  are  positive  for  anti-topo  I  auto-Abs  (based  on  Korean 
patients  enrolled  in  this  project).  We  performed  a  genome-wide  association  study 
(GWAS)  in  137  Korean  SSc  patients  and  564  sex-matched  Korean  controls  using 
the  Affymetrix  Human  SNP  Array  5.0  containing  440,734  accessible  human 
SNPs.  The  GWAS  showed  a  distinctive  peak  of  SNP  log  p-value  (p  =  7.84  x  10'15 
for  association  with  SSc)  (Figure  1).  The  peak  was  formed  with  the  SNPs 
rs31 28930,  rs7763822,  rs7764491,  rs31 17230  and  rs31 28965,  which  were 
located  in  the  region  of  HLA-DPB1  and  -DPB2  (a  psuedogene)  on  chromosome 
6p  (Figure  1,  2).  Fine  mapping  analysis  of  this  region  confirmed  that  rs31 28930, 
rs7763822  and  rs7764491  were  the  culprit  SNPs  for  associations  with  Korean 
SSc  (Figure  2,  Table  1).  Interestingly,  the  association  was  even  stronger  in 
patients  who  were  positive  for  anti-topo  I  auto-Abs  (rs31 28930,  p  =  1.70  x  1 0'22, 
OR  5.15,  95%  Cl  3.62-7.34)  (Table  1).  These  SNPs  also  were  associated  with 
the  diffuse  form  of  SSc,  but  not  the  limited  forms  of  SSc  (Supplementary  Table  1 ). 
Subtyping  of  HLA-DPB1  showed  that  HLA-DPB1*1301  (21.0%  in  SSc  vs.  5.5%  in 
controls),  DPB1*0901  (12.0%  vs.  2.6%)  and  DPB1*030101  (10.0%  vs.  4.3%) 
were  significantly  more  represented  in  anti-topo  1  positive  patients  than  controls 
(Table  2).  SNPs  corresponding  to  previous  SSc  associated  reported  genes,  such 
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as  PTPN22,  AIF1,  TNF ,  CTLA-4,  TGFB  and  CTGF,  fell  within  the  significance 
thresholds  of  10'5-10'6  advocated  for  gene-based  scans,  as  well  as  the 
Bonferroni  correction  for  multiple  comparisons18. 

To  confirm  these  results,  we  used  TaqMan  Assays  to  reexamine  the  5  SNPs 
with  the  strongest  association  from  the  Korean  GWAS  screen  in  1,106  US 
Caucasian  SSc  patients  (collected  from  US),  of  whom  16%  were  positive  for  anti- 
topo  I,  and  2,747  normal  controls,  of  which  447  were  from  our  local  collections 
(Houston,  Texas,  US.)  and  2,300  were  from  the  NIH  data  base  of  Genotype  and 
Phenotype  (dbGaP)  (http://www.ncbi.nlm.nih.gov/gap).  The  SNPs  rs7763822  and 
rs7764491  showed  highly  significant  associations  with  SSc  patients  who  were 
positive  for  anti-topo  I  auto-Abs  (p  =  7.58  x  10'17  and  4.84  x  10'16,  respectively) 
(Table  3).  The  HLA-DPB1*1301  allele  which  occurs  in  only  3%  of  US  Caucasians 
was  found  in  25%  of  anti-topo  I  positive  patients  and  conferred  the  strongest  risk  by 
exact  logistic  regression  (p=0.0001,  OR=14)  of  any  HLA  class  II  allele  (unpublished 
results).  The  SNPs  rs31 28965  and  rs31 17230  showed  strong  associations  with  SSc 
patients  who  were  positive  for  anti-centromere  auto-Abs  (p  =  3.20  x  1 0'5  and  1.12  x 
1 0'3,  respectively)  (Table  3).  In  addition,  the  pair  of  anti-topo  I  associated  SNPs  also 
showed  a  weaker  association  with  the  diffuse  form  of  SSc  (p  =  0.0070  and  0.014  for 
rs7763822  and  rs7764491  respectively)  (Supplementary  Table  2).  The  genetic 
concordance  of  the  patients  with  anti-topo  I  positivity  and  the  diffuse  form  of  SSc 
supports  clinical  observations  that  these  two  traits  within  SSc  commonly  overlap2' 3. 
However,  the  anti-centromere  associated  SNP  pairs  did  not  show  strong 
associations  with  the  limited  form  of  SSc  that  usually  occurs  in  patients  with  anti- 
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centromere  auto-Abs2;  4.  The  SNP  rs31 28930  showed  only  a  marginal  p  value  of 
0.041  in  the  limited  form  of  SSc  patients  who  were  positive  for  anti-centromere 
auto-Abs  (Supplementary  table  2).  Further  analysis  of  Caucasian  patients  who  are 
negative  for  anti-topo  I  or  anti-centromere  autoantibodies  indicated  that  they  have 
marginal  or  no  association  with  the  genotypes  of  all  the  5  SNPs  (Table3).  In 
contrast,  highly  significant  differences  were  observed  in  the  comparisons  of  patients 
with  and  without  anti-topo  I  or  anti-centromere  autoantibodies  using  corresponding 
anti-topo  I  or  anti-centromere  associated  SNP  pairs  (rs7763822  and  rs7764491 
pair,  p  <  2.86x1  O'18  or  rs31 28965  and  rs31 17230  pair,  4.77x1  O'4  respectively) 
(Supplementary  Table  3).  Interestingly,  marginally  significant  differences  (0.05  >  p  > 
0.01)  were  also  observed  in  the  comparison  of  patients  with  the  limited  and  the 
diffuse  forms  of  SSc  using  both  pairs  of  SNPs  but  the  autoantibody  associations 
were  the  strongest,  perhaps  because  HLA  alleles  are  specific  immune  response 
genes  (Table  3,  and  supplementary  Table  2). 

Finally,  we  studied  these  same  SNPs  in  African  Americans  and  Hispanics 
with  limited  numbers  of  cases  and  controls  (70  cases  vs.  90  controls  and  61 
cases  vs.  90  controls,  respectively).  Although,  the  numbers  of  subjects  examined 
in  these  two  populations  were  small,  the  SNPs  rs7764491  and  rs7763822 
showed  a  consistently  strong  association  with  anti-topo  I  positive  SSc  in  both 
African  Americans  (p=9.05x10'3,  OR=4.23,  95%  Cl=1.33-13.48  for  both  SNPs) 
and  Hispanics  (7.98x10'4,  OR=5.51,  95%  01=1.85-16.43  and  p=7.21x10'4, 
OR=5.57,  95%  01=1.87-16.62,  respectively)  (Supplementary  Table  4, 

Supplementary  Table  5). 


9 


Discussion 


For  human  complex  diseases,  such  as  SSc  and  systemic  lupus  erythematosus 
(SLE),  there  have  been  inconsistent  reports  of  genetic  associations  from  different 
study  populations  (e.g.  CTGF11: 19,  PTPN22 6;  20-22  and,  TGF-beta10;  23:  24  ).  Our 
findings  of  SSc-associated  HLA-DPB1  and  -DPB2  on  the  basis  of  autoantibodies, 
represent  the  first  replicable  report  in  Asians,  Caucasians,  African-Americans 
and  Hispanics.  These  complementary  studies  in  four  independent  populations 
provide  strong  support  for  these  identified  genetic  makers  to  confer  susceptibility 
to  subgroups  of  SSc  patients.  HLA-DPB1,  located  centromeric  to  other  HLA  class 
II  molecules,  shows  relatively  low  linkage  disequilibrium  with  other  extended 
MHC  haplotypes25.  Although  it  was  not  studied  as  extensively  as  HLA-A,  -B,  -C 
or  -DR,  it  has  similar  antigen-presentation  function  to  activate  CD4+  T  cells26.  Its 
genetic  polymorphisms  have  been  found  to  be  associated  with  chronic 
berylliosis27,  graft-vs-host  disease28,  juvenile  rheumatoid  arthrtitis29,  insulin 
dependent  diabetes  mellitus30  and  sarcoidosis31.  Some  studies  have  suggested 
the  possible  roles  of  HLA-DPB1  in  SSc,  in  the  context  of  HLA-A,  B,  C  and  DR 
molecules32-34.  However,  the  results  of  our  study  suggest  that  HLA-DPB1  may  be 
the  most  important  susceptible  gene  to  SSc,  especially  to  those  patients  with 
auto-antibodies  to  topo-1 . 
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Our  studies  indicate  that  sub-classification  and  population  origin  may  be  two 
critical  factors  in  identifying  genetic  susceptibility  markers  for  SSc.  Such  is 
evidenced  in  our  lack  of  association  in  Caucasian  SSc  with  any  of  the  5  SNPs 
found  in  the  Korean  GWAS  until  we  took  into  account  those  subgroups  defined 
by  clinically  diffuse  and  limited  forms,  but  especially  by  specific  autoantibody 
status  (anti-topo  I  and  anti-centromere  positives).  In  addition,  this  notion  may 
also  explain  previous  inconsistent  reports  of  SSc  associated  genes  in  different 
populations,  as  well  as  the  fact  that  these  genes  were  not  identified  in  our  GWAS 
of  Koreans. 

Importantly,  the  genetic  differentiations  of  SSc  with  distinctive  serological  and 
clinical  features  in  our  studies  suggest  that  SSc  should  not  be  considered  as  a 
single  disease.  While  SSc,  like  other  human  complex  diseases,  is  currently 
incurable,  sub-classification  of  SSc  on  the  basis  of  genetic  polymorphisms  for 
disease  susceptibility  may  provide  a  new  dimension  for  exploring  pathogenesis 
and  treatment  of  this  disease.  Our  results  strongly  indicate  that  sub-classification 
of  SSc  on  the  basis  of  autoantibodies  against  topo  I  and  centromeric  proteins  is 
important  in  defining  disease  susceptibility  genes  in  this  heterogeneous  disease. 


Keywords:  Systemic  sclerosis,  Genome  wide  association  study,  HLA-DPB1, 
Anti-topoisomerase  I  antibody 


11 


Acknowledgements 

This  study  was  supported  by  Korea  Health  21  R&D  Project,  Ministry  of  Health, 
Welfare  &  Family  Affairs,  R.O.K.  (A030001);  NIH/NIAMS  P50  AR054144,  NOI- 
AR02251  and  UL1RR024148,  Department  of  the  Army,  Medical  Research 
Acquisition  Activity:  PR064803  and  PR064851.  This  work  cannot  be 
accomplished  without  the  generous  provision  of  GWAS  data  and  DNAs  for 
healthy  Koreans  from  Korean  Center  for  Disease  Control.  We  also  thank  Hua 
Dong  for  her  valuable  help  in  statistical  analysis 


12 


References 


1.  Medsger,  T.A.,  Jr.  (2000).  Arthritis  and  Allied  Conditions.  In  Systemic  sclerosis 

(scleroderma):  clinical  aspects  W.J.  Koopman,  ed.  (Philadelphia,  Lippincott 
Williams  &  Wilkins),  p  1590. 

2.  Bunn,  C.C.,  and  Black,  C.M.  (1999).  Systemic  sclerosis:  an  autoantibody  mosaic.  Clin 

Exp  Immunol  117,  207-208. 

3.  Jarzabek-Chorzelska,  M.,  Blaszczyk,  M.,  Jablonska,  S.,  Chorzelski,  T.,  Kumar,  V.,  and 

Beutner,  E.H.  (1986).  Scl  70  antibody— a  specific  marker  of  systemic  sclerosis.  Br 
J  Dermatol  115,393-401. 

4.  Moroi,  Y.,  Peebles,  C.,  Fritzler,  M.J.,  Steigerwald,  J.,  and  Tan,  E.M.  (1980). 

Autoantibody  to  centromere  (kinetochore)  in  scleroderma  sera.  Proc  Natl  Acad 
SciUS  A  77,  1627-1631. 

5.  Arnett,  F.C.  (1995).  HLA  and  autoimmunity  in  scleroderma  (systemic  sclerosis).  Int 

Rev  Immunol  12,  107-128. 

6.  Gourh,  P.,  Tan,  F.K.,  Assassi,  S.,  Ahn,  C.W.,  McNeamey,  T.A.,  Fischbach,  M.,  Arnett, 

F.C.,  and  Mayes,  M.D.  (2006).  Association  of  the  PTPN22  R620W 
polymorphism  with  anti-topoisomerase  I-  and  anticentromere  antibody-positive 
systemic  sclerosis.  Arthritis  Rheum  54,  3945-3953. 

7.  Alkassab,  F.,  Gourh,  P.,  Tan,  F.K.,  McNeamey,  T.,  Fischbach,  M.,  Ahn,  C.,  Amett, 

F.C.,  and  Mayes,  M.D.  (2007).  An  allograft  inflammatory  factor  1  (AIF1)  single 
nucleotide  polymorphism  (SNP)  is  associated  with  anticentromere  antibody 
positive  systemic  sclerosis.  Rheumatology  (Oxford)  46,  1248-1251. 

8.  Tolusso,  B.,  Fabris,  M.,  Caporali,  R.,  Cuomo,  G.,  Isola,  M.,  Soldano,  F.,  Montecucco, 

C.,  Valentini,  G.,  and  Ferraccioli,  G.  (2005).  -238  and  +489  TNF-alpha  along 
with  TNF-RII  gene  polymorphisms  associate  with  the  diffuse  phenotype  in 
patients  with  Systemic  Sclerosis.  Immunol  Fett  96,  103-108. 

9.  Hudson,  F.F.,  Silver,  R.M.,  and  Pandey,  J.P.  (2004).  Ethnic  differences  in  cytotoxic  T 

lymphocyte  associated  antigen  4  genotype  associations  with  systemic  sclerosis.  J 
Rheumatol  31,  85-87. 

10.  Crilly,  A.,  Hamilton,  J.,  Clark,  C.J.,  Jardine,  A.,  and  Madhok,  R.  (2002).  Analysis  of 

transforming  growth  factor  betal  gene  polymorphisms  in  patients  with  systemic 
sclerosis.  Ann  Rheum  Dis  61,  678-681. 

11.  Fonseca,  C.,  Findahl,  G.E.,  Politicos,  M.,  Sestini,  P.,  Renzoni,  E.A.,  Holmes,  A.M., 

Spagnolo,  P.,  Pantelidis,  P.,  Feoni,  P.,  McHugh,  N.,  et  al.  (2007).  A 
polymorphism  in  the  CTGF  promoter  region  associated  with  systemic  sclerosis.  N 
Engl  J  Med  357,  1210-1220. 

12.  Chifflot,  H.,  Fautrel,  B.,  Sordet,  C.,  Chatelus,  E.,  and  Sibilia,  J.  (2008).  Incidence  and 

prevalence  of  systemic  sclerosis:  a  systematic  literature  review.  Semin  Arthritis 
Rheum  37,  223-235. 

13.  Mayes,  M.D.,  Facey,  J.V.,  Jr.,  Beebe-Dimmer,  J.,  Gillespie,  B.W.,  Cooper,  B.,  Faing, 

T.J.,  and  Schottenfeld,  D.  (2003).  Prevalence,  incidence,  survival,  and  disease 
characteristics  of  systemic  sclerosis  in  a  large  US  population.  Arthritis  Rheum  48, 
2246-2255. 


13 


14.  (1980).  Preliminary  criteria  for  the  classification  of  systemic  sclerosis  (scleroderma). 

Subcommittee  for  sclerodenna  criteria  of  the  American  Rheumatism  Association 
Diagnostic  and  Therapeutic  Criteria  Committee.  Arthritis  Rheum  23,  581-590. 

15.  Frazer,  K.A.,  Ballinger,  D.G.,  Cox,  D.R.,  Hinds,  D.A.,  Stuve,  L.L.,  Gibbs,  R.A., 

Belmont,  J.W.,  Boudreau,  A.,  Hardenbol,  P.,  Leal,  S.M.,  et  al.  (2007).  A  second 
generation  human  haplotype  map  of  over  3.1  million  SNPs.  Nature  449,  851-861. 

16.  Barrett,  J.C.,  Fry,  B.,  Mailer,  J.,  and  Daly,  M.J.  (2005).  Haploview:  analysis  and 

visualization  of  LD  and  haplotype  maps.  Bioinformatics  21, 263-265. 

17.  Purcell,  S.,  Neale,  B.,  Todd-Brown,  K.,  Thomas,  L.,  Ferreira,  M.A.,  Bender,  D., 

Mailer,  J.,  Sklar,  P.,  de  Bakker,  P.I.,  Daly,  M.J.,  et  al.  (2007).  PLINK:  a  tool  set 
for  whole-genome  association  and  population-based  linkage  analyses.  Am  J  Hum 
Genet  81,  559-575. 

18.  Thomas,  D.C.,  and  Clayton,  D.G.  (2004).  Betting  odds  and  genetic  associations.  J 

Natl  Cancer  Inst  96,  421-423. 

19.  Rueda,  B.,  Simeon,  C.,  Hesselstrand,  R.,  Herrick,  A.,  Worthington,  J.,  Ortego- 

Centeno,  N.,  Riemekasten,  G.,  Fonollosa,  V.,  Vonk,  M.C.,  van  den  Hoogen,  F.H., 
et  al.  (2008).  A  large  multicenter  analysis  of  CTGF  -945  promoter  polymorphism 
does  not  confirm  association  with  Systemic  Sclerosis  susceptibility  or  phenotype. 
Ann  Rheum  Dis. 

20.  Kyogoku,  C.,  Langefeld,  C.D.,  Ortmann,  W.A.,  Lee,  A.,  Selby,  S.,  Carlton,  V.E., 

Chang,  M.,  Ramos,  P.,  Baechler,  E.C.,  Batliwalla,  F.M.,  et  al.  (2004).  Genetic 
association  of  the  R620W  polymorphism  of  protein  tyrosine  phosphatase  PTPN22 
with  human  SLE.  Am  J  Hum  Genet  75,  504-507. 

21.  Viken,  M.K.,  Amundsen,  S.S.,  Kvien,  T.K.,  Boberg,  K.M.,  Gilboe,  I.M.,  Lilleby,  V., 

Sollid,  L.M.,  Forre,  O.T.,  Thorsby,  E.,  Smerdel,  A.,  et  al.  (2005).  Association 
analysis  of  the  1858C>T  polymorphism  in  the  PTPN22  gene  in  juvenile 
idiopathic  arthritis  and  other  autoimmune  diseases.  Genes  Immun  6,  271-273. 

22.  Wipff,  J.,  Allanore,  Y.,  Kahan,  A.,  Meyer,  O.,  Mouthon,  L.,  Guillevin,  L.,  Pierlot,  C., 

Glikmans,  E.,  Bardin,  T.,  Boileau,  C.,  et  al.  (2006).  Lack  of  association  between 
the  protein  tyrosine  phosphatase  non-receptor  22  (PTPN22)*620W  allele  and 
systemic  sclerosis  in  the  French  Caucasian  population.  Ann  Rheum  Dis  65,  1230- 
1232. 

23.  Lee,  E.B.,  Kim,  J.Y.,  Lee,  Y.J.,  Abdallah,  A.,  Lympany,  P.,  and  Song,  Y.W.  (2004). 

Transforming  growth  factor-beta  1  polymorphisms  in  Korean  patients  with 
systemic  sclerosis.  Tissue  Antigens  63,  491-495. 

24.  Sugiura,  Y.,  Banno,  S.,  Matsumoto,  Y.,  Niirni,  T.,  Yoshinouchi,  T.,  Hayami,  Y., 

Naniwa,  T.,  and  Ueda,  R.  (2003).  Transforming  growth  factor  betal  gene 
polymorphism  in  patients  with  systemic  sclerosis.  J  Rheumatol  30,  1520-1523. 

25.  Begovich,  A.B.,  McClure,  G.R.,  Suraj,  V.C.,  Helmuth,  R.C.,  Fildes,  N.,  Bugawan, 

T.L.,  Erlich,  H.A.,  and  Klitz,  W.  (1992).  Polymorphism,  recombination,  and 
linkage  disequilibrium  within  the  HLA  class  II  region.  J  Immunol  148,  249-258. 

26.  McCanlies,  E.C.,  Kreiss,  K.,  Andrew,  M.,  and  Weston,  A.  (2003).  HLA-DPB1  and 

chronic  beryllium  disease:  a  HuGE  review.  Am  J  Epidemiol  157,  388-398. 

27.  Richeldi,  L.,  Sorrentino,  R.,  and  Saltini,  C.  (1993).  HLA-DPB1  glutamate  69:  a 

genetic  marker  of  beryllium  disease.  Science  262,  242-244. 


14 


28.  Shaw,  B.E.,  Potter,  M.N.,  Mayor,  N.P.,  Pay,  A.L.,  Smith,  C.,  Goldman,  J.M., 

Prentice,  H.G.,  Marsh,  S.G.,  and  Madrigal,  J.A.  (2003).  The  degree  of  matching  at 
HLA-DPB 1  predicts  for  acute  graft-versus-host  disease  and  disease  relapse 
following  haematopoietic  stem  cell  transplantation.  Bone  Marrow  Transplant  3 1 , 
1001-1008. 

29.  Begovich,  A.B.,  Bugawan,  T.L.,  Nepom,  B.S.,  Klitz,  W.,  Nepom,  G.T.,  and  Erlich, 

H.A.  (1989).  A  specific  HLA-DP  beta  allele  is  associated  with  pauciarticular 
juvenile  rheumatoid  arthritis  but  not  adult  rheumatoid  arthritis.  Proc  Natl  Acad 
SciUS  A  86,  9489-9493. 

30.  Cruz,  T.D.,  Valdes,  A.M.,  Santiago,  A.,  Frazer  de  Llado,  T.,  Raffel,  L.J.,  Zeidler,  A., 

Rotter,  J.I.,  Erlich,  H.A.,  Rewers,  M.,  Bugawan,  T.,  et  al.  (2004).  DPB1  alleles  are 
associated  with  type  1  diabetes  susceptibility  in  multiple  ethnic  groups.  Diabetes 
53,2158-2163.  ' 

31.  Lympany,  P.A.,  Petrek,  M.,  Southcott,  A.M.,  Newman  Taylor,  A.J.,  Welsh,  K.I.,  and 

du  Bois,  R.M.  (1996).  HLA-DPB  polymorphisms:  Glu  69  association  with 
sarcoidosis.  Eur  J  Immunogenet  23,  353-359. 

32.  Kuwana,  M.,  Inoko,  H.,  Kameda,  H.,  Nojima,  T.,  Sato,  S.,  Nakamura,  K.,  Ogasawara, 

T.,  Hirakata,  M.,  Ohosone,  Y.,  Kaburaki,  J.,  et  al.  (1999).  Association  of  human 
leukocyte  antigen  class  II  genes  with  autoantibody  profiles,  but  not  with  disease 
susceptibility  in  Japanese  patients  with  systemic  sclerosis.  Intern  Med  38,  336- 
344. 

33.  Tikly,  M.,  Rands,  A.,  McHugh,  N.,  Wordsworth,  P.,  and  Welsh,  K.  (2004).  Human 

leukocyte  antigen  class  II  associations  with  systemic  sclerosis  in  South  Africans. 
Tissue  Antigens  63,  487-490. 

34.  Reveille,  J.D.,  Brady,  J.,  MacLeod-St  Clair,  M.,  and  Durban,  E.  (1992).  HLA-DPB  1 

alleles  and  autoantibody  subsets  in  systemic  lupus  erythematosus,  Sjogren's 
syndrome  and  progressive  systemic  sclerosis:  a  question  of  disease  relevance. 
Tissue  Antigens  40,  45-48. 


15 


Figures  legends 

Figure  1.  Identification  of  the  major  locus  associated  with  SSc  in  genome  wide 
scan  of  Koreans.  A  total  of  440,734  SNPs  were  evaluated  in  133  patients  with 
SSc  and  557  healthy  controls.  A)  Quantile-quantile  plot  compares  the 
distributions  of  observed  P-values  with  expected  null  P  values.  B)  Distribution  of 
-Logio  P  values  are  plotted  against  chromosomes. 

Figure  2.  HLA-DPB  as  a  culprit  region  for  susceptibility  to  systemic  sclerosis.  - 
LogioP  values  are  depicted  around  the  HLA-DPB  region.  All  genes  in  the  region 
are  also  displayed  above  the  linkage  disequilibrium  (LD)  map.  LD  between  pairs 
of  SNPs  is  depicted  with  linkage  blocks  in  which  bright  red  block  represents 
disequilibrium  coefficient  (D’)  of  1.0  and  white  block  D’  of  0.0.  Orange  bar  above 
the  LD  map  represents  high  LD  SNPs  which  show  significant  association  (r2>0.8) 
with  culprit  SNPs. 
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Figure  1 
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Figure  2 
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Table  1.  Association  of  SNPs  rs31 28930,  rs7764491  and  rs31 28965  with  Korean  SSc 
patients  with  and  without  autoantibodies  to  DNA  topoisomerase  I  and  centromere  protein. 


Alleles 

(Mninor  allele) 

Systemic  sclerosis 

Controls 

Anti-topoisomerase  1  Abb 

Anti-centromere  Ab 

Total 

(n=133) 

(n=557) 

(+) 

(n=79) 

(-) 

(n=48) 

(+) 

(n=16) 

(-) 

(n=  117) 

rs31 28930  (A) 

Frequency 

0.47 

0.15 

0.13 

0.39 

0.36 

0.15 

P-valuea 

1.70x10-22 

0.93 

0.705 

2.38x1  O'17 

7.84x1 0-15 

- 

OR  (95%  Cl)a 

5.15(3.62-7.34) 

0.97  (0.54-1.76) 

0.82(0.28-2.35) 

3.63(2.66-4.95) 

3.17(2.35-4.28) 

rs7764491  (C) 

Frequency 

0.21 

0.073 

0.063 

0.17 

0.16 

0.057 

P-valuea 

1.1 8x1 0'11 

0.51 

0.886 

2.13x10-9 

2.13x1O-0 

- 

OR  (95%  Cl)a 

4.40(2.78-6.98) 

1.31(0.58-2.95) 

1.11(0.26-4.76) 

3.44(2.25-5.26) 

3.13(2.06-4.74) 

rs31 28965  (A) 

Frequency 

0.26 

0.073 

0.063 

0.21 

0.19 

0.093 

P-valuea 

7.78x1 0-10 

0.506 

0.553 

1.44x1 0-7 

2.29x1 0-6 

- 

OR(95%  Cl)a 

3.40(2.26-5.12) 

0.76  (0.34-1.69) 

0.65(0.15-2.75) 

2.64(1.82-3.83) 

2.36(1.64-3.40) 

a. Each  subset  versus  controls  b. Missing  information  of  anti- 
topoisomerase  I  antibody  status  in  6  patients;  Ab,  antibody;  Cl, 
confidence  interval 
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Table  2  .  Association  of  HLA-DPB1  allelic  subtypes  with  systemic  sclerosis  in  Koreans. 


H  LA  alleles 

Systemic  sclerosis 

Controls 

Anti-topoisomerase-l  Abb 

Anti-centromere  Ab 

Total 

(n=1 37) 

(n=548) 

(+) 

(n=82) 

(-) 

(n=49) 

(+) 

(n=16) 

(-) 

(n=121) 

*1301 

Frequency 

0.21 

0.05 

0.031 

0.17 

0.15 

0.055 

P-valuea 

4.05x1 0-12 

0.876 

0.562 

3.25x10-9 

7.61x10-8 

- 

OR  (95%  Cl)a 

4.52(2.86-7.14) 

0.93(0.36-2.37) 

0.56(0.074-4.15) 

3.42(2.23-5.24) 

3.04(1.99-4.63) 

*0901 

Frequency 

0.12 

0.01 

0.031 

0.087 

0.08 

0.026 

P-valuea 

2.44x1  O'8 

0.325 

0.868 

7.55x1 0-6 

2.55x1 0-5 

- 

OR  (95%  Cl)a 

4.82(2.64-8.82) 

0.38(0.05-2.82) 

1.19(0.16-8.99) 

3.50(1.96-6.24) 

3.21(1.82-5.69) 

*030101 

Frequency 

0.1 

0.02 

0.031 

0.083 

0.077 

0.043 

P-valuea 

9.47x1 0-4 

0.283 

0.748 

0.01 

0.0214 

- 

OR  (95%  Cl)a 

2.58(1.44-4.61) 

0.47(0.11-1.94) 

0.72(0.096-5.39) 

2.01(1.17-3.46) 

1.85(1.09-3.16) 

a. Each  subset  versus  controls  b. Missing  information  of  anti-topoisomerase  I  antibody 


status  in  6  patients;  Ab,  antibody;  Cl,  confidence  interval 
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Table  3.  Association  of  SNPs  of  HLA -DPB1  and  -DPB2  with  systemic  sclerosis  in  Caucasians. 


Alleles 

Systemic  sclerosis15 

Controls 

(Minor  allele) 

Anti-topoisomerase  1  Ab 

Anti-centromere  Ab 

Total 

(+) 

(n=183) 

(-) 

(n=9 1 7) 

(+) 

(n=316) 

(-) 

(n=784) 

(n= 1 , 1 07) 

(n=2,731) 

rs31 28930  (A) 

Frequency 

0.31 

0.28 

0.3 

0.28 

0.29 

0.27 

P-valuea 

0.117 

0.295 

0.0564 

0.47 

0.14 

- 

OR  (95%  Cl)a 

1.20(0.95-1.51) 

1.01(0.95-1.20) 

1.19(1.00-1.43) 

1.05(0.92-1.19) 

1.09(0.97-1.21) 

- 

rs7764491  (C) 

Frequency 

0.14 

0.03 

0.025 

0.058 

0.049 

0.043 

P-valuea 

4.84x10-16 

0.0175 

0.036 

0.0155 

0.29 

OR  (95%  Cl)a 

3.56(2.57-4.94) 

0.70(0.52-0.94) 

0.58(0.35-0.97) 

1.36(1.06-1.75) 

1.14(0.90-1.44) 

- 

rs7763822  (T) 

Frequency 

0.14 

0.031 

0.024 

0.059 

0.049 

0.043 

P-valuea 

7.58x1 0-17 

0.0214 

0.0227 

0.0076 

0.236 

- 

OR  (95%  Cl)a 

3.65(2.64-5.04) 

0.71(0.53-0.95) 

0.55(0.32-0.93) 

1.40(1.09-1.80) 

1.15(0.91-1.46) 

- 

rs31 28965  (A) 

Frequency 

0.14 

0.19 

0.25 

0.16 

0.18 

0.18 

P-valuea 

0.0237 

0.394 

3.20x1  O'5 

0.01 

0.99 

- 

OR  (95%  Cl)a 

0.70(0.52-0.96) 

1.06(0.93-1.21) 

1.50(1.24-1.82) 

0.82(0.70-0.95) 

1.00  (0.88-1.14) 

- 

rs31 17230  (G) 

Frequency 

0.17 

0.26 

0.29 

0.22 

0.24 

0.24 

P-valuea 

0.00765 

0.055 

1.12x10-3 

0.32 

0.42 

- 

OR  (95%  Cl)a 

0.69(0.52-0.91) 

1.13(1.00-1.28) 

1.36(1.13-1.63) 

0.93(0.81-1.07) 

1.05(0.93-1.18) 

- 

a. Each  subset  versus  controls  b.  Missing  information  of  autoantibody  status  in  6 


patients;  Ab,  antibody;  Cl  confidence  interval 
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Supplementary  Table  1.  Associaiton  of  rs312930,  rs7764491  and  rs3128965  with 


subsets  of  systemic  sclerosis  in  Koreans. 


SNPs 

Systemic  sclerosis 

Controls 

(Minor  allele) 

Skin  subset 

Limited  Diffuse 

(n=67)  (n=66) 

Limited/Cent+ 

(n=12) 

Diffuse/To  po+ 

(n=45) 

rs31 28930  (A) 

Frequency 

0.33 

0.39 

0.13 

0.49 

0.15 

P-valuea 

1.67x10-7 

1.1 5x1 0-11 

0.742 

3.33x1  O'16 

OR  (95%  Cl)a 

2.79(1.88-4.15) 

3.59(2.44-5.29) 

0.81(0.24-2.76) 

5.46(3.50-8.51) 

rs7764491(C) 

Frequency 

0.11 

0.2 

0.083 

0.26 

0.057 

P-valuea 

0.0123 

5.28x10-10 

0.576 

1.77x1  O’12 

OR  (95%  Cl)a 

2.10(1.16-3.81) 

4.29(2.62-7.03) 

1.52(0.35-6.59) 

5.73(3.35-9.80) 

rs31 28965  (A) 

Frequency 

0.22 

0.17 

0.042 

0.22 

0.093 

P-valuea 

1.29x10-5 

3.68x1 0-3 

0.387 

1.09x1  O'4 

OR  (95%  Cl)a 

2.68(1.70-4.24) 

2.05(1.25-3.36) 

0.42(0.056-3.16) 

2.78(1.62-4.74) 

a. Each  subset  versus  controls;  Ab,  antibody;  Cent,  anti-centromere  antibody;  Cl 


confidence  interval;  Topo,  anti-topoisomerase  I  antibody 
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Supplementary  Table  2.  Association  of  SNPs  of  HLA-DPB1  and  -DPB2  with 
subsets  of  systemic  sclerosis  in  Caucasians 


SNPs 

Systemic  sclerosis 

Controls 

(Minor  allele) 

Skin  subset 

Limited  Diffuse 

(n=654)  (n=419) 

Limited/Cent+ 

(n=277) 

Diffuse/Topo+ 

(n=88) 

(n=2,731) 

rs31 28930  (A) 

Frequency 

0.28 

0.28 

0.31 

0.33 

0.27 

P-valuea 

0.238 

0.376 

0.041 

0.072 

- 

OR  (95%  Cl)a 

1.08(0.94-1.24) 

1.08(0.91-1.27) 

1.22(1.01-1.48) 

1.34(0.97-1.85) 

- 

rs7764491(C) 

Frequency 

0.04 

0.62 

0.025 

0.14 

0.043 

P-valuea 

0.647 

0.014 

0.0485 

4.76x10-10 

- 

OR  (95%  Cl)a 

0.93(0.68-1.27) 

1.47(1.08-2.02) 

0.58(0.34-1.00) 

3.73(2.39-5.81) 

- 

rs7763822(T) 

Frequency 

0.04 

0.064 

0.024 

0.15 

0.043 

P-valuea 

0.632 

0.007 

0.0299 

3.99X10-11 

- 

OR  (95%  Cl)a 

0.93(0.68-1.26) 

1.52(1.12-2.07) 

0.54(0.31-0.95) 

3.93(2.54-6.08) 

- 

rs31 28965  (A) 

Frequency 

0.2 

0.16 

0.26 

0.13 

0.18 

P-valuea 

0.22 

0.0765 

1.71x10-5 

0.047 

- 

OR  (95%  Cl)a 

1.1(0.94-1.28) 

0.84(0.69-1.02) 

1.56(1.27-1.90) 

0.63(0.40-1.00) 

- 

rs31 17230(G) 

Frequency 

0.26 

0.21 

0.3 

0.15 

0.24 

P-valuea 

0.046 

0.15 

6.9x1 0-4 

6.08x10-3 

- 

OR  (95%  Cl)a 

1.15(1.00-1.33) 

0.88(0.73-1.04) 

1.40(1.15-1.70) 

0.55(0.36-0.84) 

- 

a. Each  subset  versus  controls;  Ab,  antibody;  Cent,  anti-centromere  antibody;  Cl 


confidence  interval;  Topo,  anti-topoisomerase  I  antibody 
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Supplementary  Table  3.  Comparison  of  Caucasian  SSc  patients  with  and  without 


autoantibodies  to  topo  I  or  centromere  protein. 


Alleles 

(Minor  allele) 

Anti-topoisomerase  1  Ab 

Anti-centromere  Ab 

p-valuesb 

OR(95%  Cl) 

(+) 

(n=183) 

(-) 

(n=917) 

p-valuesa 

OR(95%  Cl) 

(+) 

(n=316) 

(-) 

(n=784) 

rs31 28930  (A) 

Frequency 

0.31 

0.28 

0.33 

0.30 

0.28 

0.212 

1.13(0.88-1.44) 

1.14(0.93-1.39) 

rs7764491  (C) 

Frequency 

0.14 

0.030 

2.86x10'18 

0.025 

0.058 

1 .5x1  O'3 

5.11(3.42-7.63) 

0.43(0.25-0.73) 

rs7763822  (T) 

Frequency 

0.14 

0.031 

8.61x10'19 

0.024 

0.059 

5.60x1 0-4 

5.12(3.46-7.68) 

0.39(0.22-0.68) 

rs31 28965  (A) 

Frequency 

0.14 

0.19 

0.0116 

0.25 

0.16 

1.07x10-7 

0.66(0.48-0.91) 

1.84(1.46-2.30) 

rs31 17230  (G) 

Frequency 

0.17 

0.26 

7.35x1 0-4 

0.29 

0.22 

4.77x1 0-4 

0.61(0.45-0.81) 

1.45(1.18-1.79) 

a.  Comparison  between  anti-topoisomerase  I  antibody  +  and  b.  Comparison  between  anti¬ 


centromere  antibody  +  and 
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Supplementary  Table  4.  Association  of  the  SNPs  with  SSc  patients  with  and  without 
autoantibodies  to  DNA  topoisomerase  I  and  centromere  protein  in  African  Americans  . 


SNPs 

Systemic  sclerosis 

Controls 

(Minor  allele) 

Anti-topisomerase-l  Ab 

Anti-centromere  Ab 

Total 

(+) 

(-) 

(+) 

(-) 

(n=9) 

(n=61) 

(n=5) 

(n=65) 

(n=70) 

(n=90) 

rs31 28930  (A) 

Frequency 

0.5 

0.55 

0.3 

0.56 

0.54 

0.44 

P-valuea 

0.615 

0.0588 

0.391 

0.0325 

0.0637 

- 

OR  (95%  Cl)a 

1.28(0.49-3.38) 

1.56(0.98-2.48) 

0.55(0.14-2.19) 

1.64(1.04-2.59) 

1.52(0.98-2.38) 

- 

rs7764491(C) 

Frequency 

0.28 

0.14 

0.1 

0.16 

0.16 

0.083 

P-valuea 

9.05x10-3 

0.109 

0.853 

0.0297 

0.036 

- 

OR  (95%  Cl)a 

rs7763822(T) 

4.23(1.33-13.48) 

1.82(0.87-3.79) 

1.22(0.14-10.31) 

2.16(1.07-4.37) 

2.09(1.04-4.19) 

- 

Frequency 

0.28 

0.14 

0.1 

0.16 

0.16 

0.083 

P-valuea 

9.05X10-3 

0.121 

0.853 

0.0339 

0.0405 

- 

OR  (95%  Cl)a 

rs31 28965  (A) 

4.23(1.33-13.48) 

1.78(0.85-3.72) 

1.22(0.14-10.31) 

2.12(1.05-4.29) 

2.05(1.02-4.12) 

- 

Frequency 

0 

0.067 

0 

0.063 

0.058 

0.068 

P-valuea 

0.253 

0.959 

0.393 

0.844 

0.713 

- 

OR  (95%  Cl)a 

rs31 17230(G) 

- 

0.98(0.39-2.47) 

- 

0.91(0.36-2.30) 

0.84(0.33-2.12) 

- 

Frequency 

0.28 

0.41 

0.2 

0.41 

0.39 

0.38 

P-valuea 

0.383 

0.628 

0.247 

0.649 

0.844 

- 

OR  (95%  Cl)a 

0.62(0.21-1.82) 

1.12(0.70-1.80) 

0.40(0.83-1.96) 

1.11(0.70-1.77) 

1.05(0.66-1.65) 

- 

a. Each  subset  versus  controls;  Ab,  antibody;  Cl  confidence  interval 


GWAS  of  systemic  sclerosis 


Zhou  Page  26 


Supplementary  Table  5.  Association  of  the  SNPs  with  SSc  patients  with  and  without 
autoantibodies  to  DNA  topoisomerase  I  and  centromere  protein  in  Hispanics. 


SNPs 

Systemic  sclerosis 

Controls 

(Minor  allele) 

Anti-topisomerase-l  Ab 

Anti-centromere  Ab 

Total 

(+) 

(n=1 7) 

(-) 

(n=44) 

(+) 

(n=12) 

(-) 

(n=49) 

(n=61) 

(n=90) 

rs31 28930(A) 

Frequency 

0.29 

0.38 

0.29 

0.37 

0.35 

0.24 

P-valuea 

0.493 

0.02 

0.571 

0.0236 

0.0324 

- 

OR(  95%  Cl)a 

1.33(0.59-3.00) 

1.91(1.10-3.33) 

1.31(0.51-3.38) 

1.85(1.01-3.17) 

1.74(1.05-2.89) 

- 

rs7764491(C) 

Frequency 

0.21 

0.034 

0 

0.1 

0.082 

0.045 

P-valuea 

7.98x1 0-4 

0.676 

0.289 

0.066 

0.185 

- 

OR  (95%  Cl)a 

5.51(1.85-16.43) 

0.75(0.19-2.90) 

- 

2.42(0.92-6.34) 

1.90(0.73-4.95) 

- 

rs7763822(T) 

Frequency 

0.21 

0.034 

0 

0.01 

0.082 

0.044 

P-valuea 

7.21x10-4 

0.688 

0.292 

0.0623 

0.177 

- 

OR  (95%  Cl)a 

5.57(1.87-16.62) 

0.76(0.20-2.93) 

- 

2.44(0.93-6.41) 

1.92(0.74-5.01) 

- 

rs31 28965(A) 

Frequency 

0.029 

0.2 

0.25 

0.13 

0.16 

0.16 

P-valuea 

0.0429 

0.381 

0.278 

0.532 

0.904 

- 

OR  (95%  Cl)a 

0.16(0.020-1.20) 

1.34(0.70-2.59) 

1.73(0.63-4.77) 

0.80(0.39-1.62) 

0.96(0.51-1.82) 

- 

rs31 17230(G) 

Frequency 

0.088 

0.35 

0.25 

0.29 

0.28 

0.21 

P-valuea 

0.0972 

0.013 

0.656 

0.159 

0.173 

- 

OR  (95%  Cl)a 

0.36(0.11-1.26) 

2.04(1.16-3.61) 

1.25(0.46-3.38) 

1.50(0.85-2.65) 

1.45(0.85-2.48) 

- 

a. Each  subset  versus  controls;  Ab,  antibody;  Cl  confidence  interval 
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ABSTRACT 

Objective:  SPARC  (secreted  protein,  acidic  and  rich  in  cysteine)  is  a  matri cellular  protein, 
which  along  with  other  extracellular  matrix  (ECM  )  components  including  collagens,  is  commonly 
over-expressed  in  fibrotic  diseases.  The  purpose  of  this  study  was  to  examine  whether  inhibition  of 
SPARC  can  regulate  collage  expression  in  vitro  and  in  vivo,  and  subsequently  attenuate  fibrotic 
stimulation  by  bleomycin  in  mouse  skin  and  lungs. 

M  ethods:  In  in  vitro  studies,  skin  fibroblasts  obtained  from  aTgfbrl  knock-in  mouse  (TBR1CA; 
Cre-ER)  were  transfected  with  Sparc  siRNA.  Gene  and  protein  expressions  of  the  C ol  1A 2  and  the 
Ctgf  were  examined  by  real-time  RT-PCR  and  western  bloting,  respectively.  In  in  vivo  studies, 
C57BL/6  mice  were  induced  for  skin  and  lung  fibrosis  by  bleomycin  and  followed  by  Sparc  siRNA 
treatment  through  subcutaneous  injection  and  intratracheal  instillation,  respectively.  The 
pathological  changes  of  skin  and  lungs  were  assessed  by  hematoxylin  and  eosin  (HE)  and  Masson's 
trichrome  stains.  The  expression  changes  of  collagen  in  the  tissues  were  assessed  by  the  real-time 
RT-PCR  and  non-crosslinked  fibrillar  collagen  content  assays. 

Results:  Sparc  siRNA  significantly  reduced  gene  and  protein  expression  of  collagen  type  1  in 
fibroblasts  obtained  from  the  TBR1CA;  Cre-ER  mouse  that  was  induced  for  constitutive! y  active 
TGF-p  receptor  I.  Skin  and  lung  fibrosis  induced  by  bleomycin  was  markedly  reduced  by  treatment 
with  Sparc  siRNA.  The  anti-fibrotic  effect  of  Sparc  siRNA  in  vivo  was  accompanied  by  an 
inhibition  of  Ctgf  expression  in  these  same  tissues. 

Conclusion:  Specific  inhibition  of  Sparc  effectively  reduced  fibrotic  changes  in  vitro  and  in 
vivo.  SPARC  inhibition  may  represent  a  potential  therapeutic  approach  to  fibrotic  diseases. 
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INTRODUCTION 

Fibrosis  is  a  general  pathological  process  in  which  excessive  deposition  of  extracellular 
matrix  (ECM  )  occurs  in  the  tissues.  It  is  currently  untreatable.  Although  therapeutic  uses  of  some 
anti-inflammatory  and  immunosuppressive  agents  such  as  colchicines,  interferon-gamma, 
corticosteroids  and  cyclophosphamide  have  been  reported,  many  of  these  approaches  have  not 
proven  successful  (1-3).  Recently,  SPARC  (secreted  protein,  acidic  and  rich  in  cysteine),  a 
matricellular  component  of  the  ECM,  has  been  reported  as  a  bio-marker  for  fibrosis  in  multiple 
fibrotic  diseases,  such  as  interstitial  pulmonary  fibrosis,  renal  interstitial  fibrosis,  cirrhosis, 
atherosclerotic  lesions  and  scleroderma  or  systemic  sclerosis  (SSc)  (4-9).  Notably,  increased 
expression  of  SPARC  has  been  observed  in  affected  skin  and  circulation  of  patients  with  SSc 
(10,11),  a  devastating  disease  of  systemic  fibrosis,  as  well  as  in  cultured  dermal  fibroblasts 
obtained  from  SSc  skin  (8,9). 

SPARC,  also  called  osteonectin  or  BM-40,  is  an  important  mediator  of  cell-matrix 
interaction  (12).  Increasing  evidence  indicates  that  SPARC  may  play  an  important  role  in  tissue 
fibrosis.  In  addition  to  its  higher  expression  level  in  the  tissues  of  fibrotic  diseases,  SPARC  has 
shown  a  capacity  to  stimulate  the  transforming  growth  factor  beta  (TGF-p)  signaling  system  (13). 
Inhibition  of  SPARC  attenuates  the  profibrotic  effect  of  exogenous  TGF-p  in  cultured  human 
fibroblasts  (14).  Moreover,  in  animal  studies,  SPARC-null  mice  display  a  diminished  amount  of 
pulmonary  fibrosis  compared  with  control  mice  after  exposure  to  bleomycin,  a  chemotherapeutic 
antibiotic  with  a  profibrotic  effect  (15).  These  observations  suggest  that  SPARC  is  a  potential 
bio-target  for  anti-fibrotic  therapy.  The  studies  described  here  aimed  to  explore  the  feasibility  of 
inhibition  of  SPARC  with  siRNA  to  counter  fibrotic  processes  both  in  vitro  and  in  vivo  using 
murine  models. 

MATERIALS  AND  METHODS 


Fibroblast  cell  lines  from  Tgfbrl  knock-in  mouse.  Constitutively  activated  T gfbrl  mice, 
which  recapitulated  clinical,  histological,  and  biochemical  features  of  human  SSc,  have  been 
reported  previously  (16).  They  are  termed  TBR1CA;  Cre-ER  mice  and  harbor  both  the  DNA  for  an 
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inducible  constitutively  activeTGFp  receptor  I  (T G FpR I )  mutation  targeted  to  the  ROSA  locus,  and 
a  Cre-ER  transgene  driven  by  a  Coll  fibroblast-specific  promoter.  Administration  of 
4-hydroxytamoxifen  (4-OHT)  2  weeks  after  birth  activates  the  expression  of  constitutively  active 
Tgfbrl  in  the  mice  (16).  Fibroblasts  were  derived  from  skin  biopsy  specimens  of  these  mice.  The 
cultures  were  maintained  in  DM  EM  with  10%  FCS  and  supplemented  with  antibiotics  (50  U/ml 
penicillin  and  50  pg/ml  streptomycin).  Fifth-passage  fibroblast  cells  were  seeded  at  a  density  of  5  x 
105  cells  in  25-cm2  flasks  and  grown  until  confluence.  Experiments  were  performed  in  triplicates. 

Transient  transfection  with  siRNA  in  fibroblasts.  Double-stranded  ON -TARGET  plus 
siRNAs  of  murine  Sparc  and  Ctgf  were  purchased  from  Dharmacon,  Inc.  (Lafayette,  CO).  The 
corresponding  target  sequences  are  5'-  GCACCACACGUUUCUUUG  -3'  for  Sparc  and  5'- 
GCACCAGUGUGAAGACAUA  -3'  for  Ctgf,  respectively.  The  culture  medium  in  each  culture 
flask  with  confluent  fibroblasts  was  replaced  with  Opti-M  EM  I  medium  (Invitrogen,  Carlsbad,  CA) 
without  FCS  and  antibiotics.  The  fibroblasts  were  incubated  for  24  hours  and  transfected  with  Sparc 
siRNA  or  Ctgf  siRNA  in  a  concentration  of  100  nmol/L,  using  DharmaFFCF™  1  siRNA 
Transfection  Reagent  (Dharmacon).  Fibroblasts  with  Non-Targeting  siRNA  (Dharmacon)  treatment 
were  used  as  negative  controls.  After  24  hours,  the  culture  medium  was  replaced  with  DM  EM  .  The 
cells  transfected  with  siRNA  were  examined  after  72  hours  of  transfection  and  used  for  RNA  and 
protein  expression  analysis.  The  experiments  were  performed  in  triplicates. 

Animal  models  of  fibrosis.  C57BL/6  mice  of  about  20  grams  were  purchased  from 
J  ackson  Laboratory  (Bar  Harbor,  M  aine).  Bleomycin  from  Teva  Parenteral  M  edicines  Inc.  (Irvine, 
CA)  was  dissolved  in  saline  and  used  in  the  mice  at  a  concentration  of  3.5  unit/kg.  Pulmonary 
fibrosis  was  induced  in  these  mice  with  once  intratracheal  instillation  of  bleomycin  in  60  pi.  For 
dermal  fibrosis,  female  C57BL/6  mice  at  6  weeks  (weighing  about  20  g)  were  treated  daily  for  4 
weeks  with  local  subcutaneous  injection  of  100  pi  bleomycin  in  the  shaved  lower  back.  Four  mice 
were  used  in  each  group.  The  animal  protocols  were  approved  by  the  Center  for  Laboratory  Animal 
Medicine  and  Care  in  the  University  of  Texas  Health  Science  Center  at  Houston  and  the 
Institutional  Animal  Use  and  Care  Committee  of  M  .D.  Anderson  CancerCenter. 

Administration  of  siRNAs  in  vivo.  For  pulmonary  fibrosis,  3  pg  of  siRNA  for  in  vivo  use 
(siSTABLE,  Dharmacon),  mixed  with  DharmaFFCF™  1  siRNA  Transfection  Reagent,  was 
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administrated  intratracheal  I  y  in  60  pi  on  days  2,  5,  12  after  bleomycin  treatment.  All  the  mice  were 
sacrificed  on  day  23  after  anesthesia,  and  the  lung  samples  were  collected.  The  left  lungs  were  fixed 
by  4%  formalin  and  used  for  further  histological  analysis.  The  right  lungs  were  minced  to  small 
pieces  and  divided  into  2  parts,  one  for  RNA  extraction  and  one  for  collagen  content  analysis.  For 
dermal  fibrosis,  the  above  siRNAs  were  injected  into  the  same  area  as  that  of  bleomycin  3  hours 
after  bleomycin  treatment  and  continued  for  4  weeks.  The  mice  were  sacrificed  on  day  29  and  the 
skin  samples  were  collected.  Saline  was  used  as  negative  control  in  both  fibrosis  studies. 

Determination  of  gene  expression  by  quantitative  RT-PCR.  Total  RNA  from  each  cell 
line  was  extracted  from  the  cultured  fibroblasts  using  RNeasy  M  ini  Kit  (Qiagen,  Valencia,  CA).  For 
mice  lung  and  skin  tissues,  the  minced  samples  were  homogenized  in  lyses  solution  (Sigma-Aldrich, 
St.  Louis,  M  0)  with  a  blender.  Then  total  RNA  was  extracted  using  Gen  Elute™  M  ammalian  Total 
RNA  M  ini  prep  Kit  (Sigma-Aldrich).  Complementary  DNA  (cDNA)  was  synthesized  using 
MultiScribe™  Reverse  Transcriptase  (Applied  Biosystems).  Quantitative  real-time  RT-PCR  was 
performed  using  an  ABI  7900  Sequence  Detector  System  (Applied  Biosystems,  Foster  City,  CA). 
The  specific  primers  and  probes  for  each  gene  ( Colla2 ,  Col3Al,  Ctgf,  and  Sparc )  were  purchased 
from  the  A ssays-on- Demand  product  line  (Applied  Biosystems).  Synthesized  cDNAs  were  mixed 
with  primers/probes  in  2  x  TaqMan  universal  PCR  buffer  and  then  assayed  on  an  ABI  7900 
sequence  detector.  The  data  obtained  from  the  assays  were  analyzed  with  SDS  2.2  software 
(Applied  Biosystems).  The  expression  level  of  each  gene  in  each  sample  was  normalized  with 
Gapdh  transcript  level. 

Western  blot  analysis.  The  cellular  lysates  extracted  from  the  cultured  fibroblasts  were 
used  for  protein  assays.  The  protein  concentration  was  determined  by  a  spectrophotometer  using 
Bradford  protein  assay  kit  (Bio-Rad  Laboratories,  Hercules,  CA).  Equal  amounts  of  protein  from 
each  sample  were  subjected  to  sodium  dodecyl  sulfate-polyacrylamide  gel  electrophoresis. 
Resolved  proteins  were  transferred  onto  PVDF  membranes  and  incubated  with  respective  primary 
antibodies,  including  anti-type  I  collagen  antibody  (Biodesign  International,  Saco,  ME),  anti-Ctgf 
antibody  (GeneTex  Inc,  San  Antonio,  TX),  and  anti-SPARC  antibody  (R&D  Systems  Inc, 
M  inneapolis,  M  N).  M  ouse  p-actin  (Alexis  Biochemicals,  San  Diego,  CA)  was  used  as  an  internal 
control.  The  secondary  antibody  was  peroxidase-conjugated  anti-rabbit,  anti-goat,  or  anti-mouse 


John  Wiley  &  Sons 


Arthritis  &  Rheumatism 


Page  7  of  20 


IgG.  Specific  proteins  were  detected  by  chemiluminescence  using  an  enhanced  chemiluminescence 
system  (Amersham,  Piscataway,  NJ).  The  intensity  of  the  bands  was  quantified  using  ImageQuant 
software  (M  olecular  Dynamics,  Sunnyvale,  CA). 

Determination  of  collagen  content.  Non-crossl inked  fibrillar  collagen  in  lung  samples 
and  skin  samples  was  measured  using  theSircol  colorimetric  assay  (Biocolor,  Belfast,  UK).  M  inced 
tissues  were  homogenized  in  0.5M  acetic  acid  with  about  1:10  ratio  of  pepsin  (Sigma-Aldrich). 

Tissues  were  weighted,  and  then  incubated  overnight  at  4°C  with  vigorous  stirring.  Digested 

samples  were  centrifuged  and  the  supernatant  was  used  for  the  analysis  with  the  Sircol  dye  reagent. 
The  protein  concentration  was  determined  using  Bradford  protein  assay  kits  and  the  collagen 
content  of  each  sample  was  normalized  to  total  protein. 

Histological  analysis.  Both  the  tissue  samples  of  lung  and  skin  were  fixed  in  4%  formalin 
and  embedded  in  paraffin.  Sections  of  5  pm  were  stained  either  with  hematoxylin  and  eosin  (HE) 
and  M  asson's  trichrome. 

Statistical  Analysis.  Results  were  expressed  as  mean  ±  SD.  The  difference  between 
different  conditions  or  treatments  was  assessed  by  Students' t- test.  A  p-value  of  less  than  0.05  was 
considered  statistically  significant. 

RESULTS 


Gene  and  protein  expression  of  Colla2,  Ctgf  and  Sparc  in  the  fibroblasts  from  TBR1CA; 
Cre-ER  mice  with  and  without  transfection  of  siRNAs 

As  measured  by  quantitative  real-time  RT-PCR,  Colla2,  Ctgf  and  Sparc  showed  increased 
expression  in  the  fibroblasts  from  TBR1CA;  Cre-ER  mice  injected  with  4-OHT,  in  which  Tgfbrl 
was  constitutive! y  active,  compared  with  those  in  the  cells  from  TBR1CA;  Cre-ER  mice  injected 
with  oil  (Figure  1).  The  fold-changes  of  each  gene  in  4-OHT-injected  TBR1CA;  Cre-ER  mice 
fibroblasts  were  3.06  ±  1.42  for  Colla2  (P  =  0.050),  4.15  ±  1.18  for  Ctgf[ P  =  0.049),  and  2.49  ± 
0.63  for  Sparc  (P  =  0.017),  respectively.  To  study  whether  inhibition  of  Sparc  induced  a  reduction 
of  collagen  in  the  fibroblasts  from  constitutive! y  active  Tgfbrl  mice,  we  transfected  Sparc  siRNA 
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into  cultured  fibroblasts  obtained  from  TBR1CA;  Cre-ER  mice  injected  with  4-OHT.  Ctgf  is  a 
down-stream  gene  in  the  TGF-p  pathway  (17-20).  Inhibition  of  ctgf  reduced  expression  of  the 
fibrotic  effect  of  TGF-p  (21).  We  used  Ctgf  siRNA  as  a  positive  control  for  inhibition  of  Ctgf  and 
collagen  expression.  Transfection  efficiency  of  siRNAs  into  fibroblasts  was  measured  using 
fluorescent  RNA  duplex  s\GLO  Green  transfection  indicator  (Dharmacon)  and  was  determined  to 
be  over  80%.  The  gene  expression  levels  from  the  Non-Targeting  siRNA  treated  fibroblasts  were 
compared  with  those  from  saline-treatment  fibroblasts,  and  no  significant  differences  were  found 
(1.05  ±  0.18-folds  for  Colla2,  1.14  ±  0.16-folds  for  Ctgf,  and  1.12  ±  0.12-folds  for  Sparc). 
Therefore,  in  the  following  in  vitro  study,  fibroblasts  with  Non-Targeting  siRNA  treatment  were 
used  as  negative  controls.  Seventy-two  hours  after  Sparc  siRNA  or  Ctgf  siRNA  transfection, 
significant  reductions  of  Sparc  (95%)  by  Sparc  siRNA  and  Ctgf  (64%)  by  Ctgf  siRNA  were 
observed  in  the  fibroblasts  (Figure  2A).  In  parallel,  Colla2  showed  decreased  expression  in  both 
siRNA  transfected  fibroblasts  (27%  and  29%  decrease  with  p  <  0.05  for  Ctgf  siRNA  and  Sparc 
siRNA,  respectively)  (Figure  1).  Western  blot  analysis  showed  a  similar  level  of  protein  reduction 
of  type  I  collagen  by  either  Sparc  siRNA  or  Ctgf  siRNA  treatment.  As  illustrated  in  Figure  2B  and 
2C,  both  Sparc  siRNA  and  Ctgf  siRNA  showed  significant  attenuation  of  collagen  type  I  in  the 
fibroblasts  (p  =  0.009  or  0.015,  respectively).  CTGF  and  SPARC  protein  levels  also  were  reduced 
by  Sparc  siRNA  (p  =  0.056  or  0.0004,  respectively)  or  Ctgf  siRNA  treatment  (p  =0.002  or  0.355). 

siRNAs  ameliorated  fibrosis  in  skin  and  reduced  inflammation  in  lungs  induced  by  bleomycin 

FIE  stains  of  mouse  skin  tissues  (Figure  3)  showed  that  4-week  injections  of  bleomycin  induced 
significant  fibrosis  in  skin  where  the  fat  cells  were  replaced  by  fiber  bundles  (Figure  3.  B), 
compared  with  normal  skin  injected  with  saline  only  (Figure  3.).  Bleomycin- injected  skin  treated 
with  Ctgf  siRNA  or  Sparc  siRNA  showed  that  most  of  the  fat  cells  still  existed  in  the  dermis 
without  prominent  fiber  bundles  (Figure  3.).  Masson's  trichrome  staining  of  the  samples  also 
showed  the  same  results.  Notably,  increased  hair  follicles  were  inconsistently  seen  in  Ctgf  siRNA- 
and  Sparc  si  R  N  A -treated  bleomycin-induced  skins. 

FIE  stain  of  mouse  lung  tissues  (Figure  4)  showed  a  significant  disruption  of  the  alveolar  units  and 
infiltration  of  inflammatory  cells  in  the  lungs  induced  by  bleomycin  (Figure  4).  Flowever,  after 
treatment  with  Ctgf  siRNA  or  Sparc  siRNA,  the  disruption  of  the  alveoli  was  improved  with  less 
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infiltrating  inflammatory  cells  (Figure  4.). 

siRNAs  attenuated  over-expression  of  collagen  and  other  fibrotic  ECM  genes  induced  by 
bleomycin  in  skin  and  lung  tissues 

Bleomycin  injection  induced  an  up-regulation  of  the  Colla2,  Col3al,  Ctgf  and  Sparc  in  both 
skin  (Figure  5A,  p  =  0.028,  0.016,  0.049  and  0.0005,  respectively)  and  lung  tissues  (Figure  5B,  p  = 
0.015,  0.005,  0.041  and  0.056,  respectively)  of  the  mice.  However,  in  Ctgf  siRNA  or  Sparc  siRNA 
treated  mice  skin  that  also  received  bleomycin  injection,  the  expression  of  the  Colla2  and  Col3al 
appeared  to  be  normal  in  skin  tissues  (Figure  5A,  p  =  0.025  and  0.003  for  each  gene  in  Ctgf  siRNA 
treatment,  and  p  =  0.031  and  0.010  in  Sparc  siRNA  treatment),  and  were  significantly  improved  in 
lung  tissues  (about  2.7-fold  reduction  for  Colla2  and  1.9-fold  reduction  for  Col3al,  compared  to 
bleomyin-injected  mice  without  siRNA  treatment,  p  <  0.05  for  both)  (Figure  5  B).  In  addition  to 
collagen  gene  expression,  the  Cr^and  the  Sparc  expression  were  significantly  reduced  by  Sparc 
siRNA  and  Ctgf  siRNA  treatment,  respectively  (Figure  5B).  In  detail,  compared  to 
bleomycin-induced  skin  and  lungs,  Sparc  siRNA  normalized  Ctgf  expression  in  both  skin  and  lungs 
(2.6-fold  reduction  in  both  with  p  =  0.1  and  0.039,  respectively).  Similarly,  Ctgf  siRNA  also 
reduced  Sparc  expression  in  skin  and  lungs  (2.9-fold  and  1.5-fold  reduction  with  p  =  0.044  and 
0.188,  respectively). 

siRNAs  reduced  the  collagen  contents  in  bleomycin-induced  mouse  skin  and  lung  samples:  T  O 
further  evaluate  anti-fibrotic  effects  of  siRNAs  on  the  fibrogenesis  of  skin  and  lung,  the  collagen 
content,  a  key  marker  in  assessing  collagen  deposition,  was  measured  in  the  collected  dermal  and 
pulmonary  samples.  Quantification  of  total  collagen  in  skin  samples  with  the  Sircol  assay  showed  a 
2.2-fold  increase  in  bleomycin-induced  skin  compared  with  saline- injected  skin  (p  =  0.0498).  Ctgf 
siRNA  treatment  reduced  the  collagen  content  to  47.6%  (p  =  0.03)  of  that  in  bleomycin-induced 
skin,  and  Sparc  siRNA  treatment  reduced  the  collagen  content  to  64.6%  (p  =  0.08)  (Figure 
6A).There  was  no  significant  difference  of  collagen  reduction  (p  =  0.08)  between  Sparc  siRNA 
treatment  and  Ctgf  siRNA  treatment. 

The  siRNA  treatments  also  showed  a  reduction  of  collagen  in  the  lung  tissues  of 
bleomycin-induced  mice  (Figure  6B).  In  bleomycin-induced  mice,  collagen  content  of  lung  tissues 
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was  3.6-fold  higher  than  that  in  control  mice  (saline-injected  mice,  p  =  0.0135).  In  Ctgf  siRNA  or 
Sparc  siRNA  treated  mice  that  also  were  bleomycin-induced,  collagen  content  of  lung  tissues  was 
68%  (p  =  0.128)  or  58%  (p  =  0.019)  of  that  in  bleomycin-induced  mice  without  siRNA  treatment. 
No  significant  difference  of  collagen  content  was  found  between  Sparc  siRNA  treatment  and  Ctgf 
siRNA  treatment  in  bleomycin-injured  lungs  (p  =  0.28). 


DISCUSSION 

Although  fibrosis  is  usually  an  irreversible  pathological  condition,  targeting  underlying  molecular 
effectors  may  reverse  an  active  status  of  the  fibrotic  process,  and  subsequently  inhibit  fibrosis.  The 
TGF-p  signaling  pathway  is  associated  with  active  fibrosis  (17).  It  begins  with  the  binding  of  the 
TGF-p  ligand  to  the  TGF-p  type  II  receptor,  which  catalyses  the  phosphorylation  of  the  type  I 
receptor  on  the  cell  membrane.  The  type  I  receptor  then  induces  the  phosphorylation  of 
receptor- regulated  SMADs  (R-SMADs)  that  bind  the  coSMAD.  The  phosphorylated 
R-SM  AD/coSM  AD  complex  enters  the  nucleus  acting  as  transcription  factors  to  regulate  target 
gene  expression  (18).  CTGF  (connective  tissue  growth  factor)  is  a  down-stream  gene  that  can  be 
activated  by  the  TGF-p  signaling  pathway  (19).  Activation  of  CTGF  is  associated  with  potent  and 
persistent  fibrotic  changes  in  the  tissues,  which  is  typically  represented  as  accumulation  of  the  EC  M 
components  including  collagens  (20). 

The  studies  described  here  utilized  the  fibroblasts  obtained  from  the  TBR1CA;  Cre-ER  mice  that 
were  induced  for  constitutively  active  TGF-p  receptor  I.  After  transfection  of  Sparc  siRNA,  the 
fibroblasts  showed  a  decreased  expression  of  Colla2  that  was  originally  over-expressed  in  the 
TBR1ca;  Cre-  ER  mice  (Figure  2).  This  phenomenon  suggests  that  SPARC  inhibition  may  interrupt 
fibrotic  TGF-p  signaling.  Although  the  specific  mechanism  for  this  suppression  is  unclear,  multiple 
previous  studies  have  demonstrated  a  mutual  regulatory  relationship  between  SPARC  and  TGF-p 
signaling  (14,21-23).  This  notion  also  is  supported  by  the  observation  of  an  over-expression  of 
Sparc  in  the  fibroblasts  of  the  BR1CA;  Cre-ER  mice  (Figure  1).  It  should  be  noted  that  the  Ctgf 
expression  in  the  fibroblasts  was  not  reduced  upon  Sparc  inhibition.  These  results  appear  to 
contradict  our  previous  report  of  parallel  inhibition  of  SPARC  and  CTGF  expression  in  human 
fibroblasts  by  SPARC  siRNA  (14).  A  possible  explanation  is  that  over-expressed  Ctgf  from 
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constitutively  activated  TGF-p  signaling  in  these  fibroblasts  may  confer  resistance  to  a 
down- regulatory  effect  from  Sparc  siRNA.  However,  such  resistance  appeared  to  have  limited 
influence  on  any  down- regulatory  effect  of  Sparc  siRNA  on  collagen  type  1,  which  suggests  that 
CTGF  is  nota  sole contri butor  to TG F-p  signaling-associated  fibrosis. 

Bleomycin  induced  fibrosis  in  mice  usually  occurs  after  inflammation  in  which  TGF-p  is 
up-regulated  (24).  Our  in  vivo  application  of  Sparc  siRNA  demonstrated  that  inhibition  of  Sparc 
significantly  reduced  inflammation  and  fibrosis  in  skin  and  lungs  induced  by  bleomycin.  In  the 
treatment  of  skin  fibrosis,  Sparc  siRNAs  reduced  fiber  bundles  accumulated  in  the  dermis  with  less 
mononuclear  cell  infiltrates  (Figure  3).  In  addition  to  histological  changes,  the  thickness  of 
bleomycin-induced  skin  treated  with  Sparc  siRNA  showed  over  50%  reduction  compared  to  that 
without  Sparc  siRNA  treatment  (data  not  shown).  The  changes  of  tissue  fibrotic  level  further  were 
confirmed  with  significantly  decreased  collagen  gene  expression  (Figure  5A).  Non-crosslinked 
fibrillar  collagen  in  the  skin  tissues  also  showed  an  average  of  35.4%  reduction  after  Sparc  siRNA 
treatment  (Figure  6A). 


In  the  treatment  of  lungs,  Sparc  siRNA  reduced  the  disruption  of  the  alveoli  induced  by  bleomycin 
(Figure  4),  which  was  accompanied  with  attenuated  gene  expression  and  protein  content  of 
collagens  as  compared  to  that  without  siRNA  treatment  (Figure  5B  and  6B).  However,  according  to 
histological  examinations  and  the  changed  levels  of  collagen  genes  and  proteins,  such  improvement 
appeared  to  be  less  perfect  as  compared  to  that  in  normal  control  mice.  On  the  other  hand,  much 
higher  levels  of  gene  expression  of  CollA2  and  Col3Al,  and  protein  content  of  collagen  were 
observed  in  bleomycin-induced  lung  tissues  when  they  were  compared  to  that  in  skin  tissues 
(5.2-fold,  6.7-fold  and  3.6-fold  increase  vs.  3.8-fold,  2.8-fold  and  2.2-fold  increase,  respectively), 
which  suggested  that  tissue  damage  and  fibrosis  in  lung  might  be  more  severe  than  that  in  skin.  In 
this  case,  treatment  of  bleomycin-induced  lung  damage  might  present  a  bigger  challenge  than  that 
of  skin,  and  the  siRNA  treatment  through  intratracheal  instillation  may  be  in  need  of  further 
optimization.  These  notions  were  supported  by  similar  findings  in  the  treatment  with  the  Ctgf 
siRNA,  a  positive  control  for  anti-fibrotic  effects. 
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Nevertheless,  Sparc  inhibition  showed  a  clear  anti-fibrotic  effect  in  bleomycin-induced  skin  and 
lung  tissues.  Notably,  these  changes  were  accompanied  with  a  significant  down  regulation  of  Ctgf 
that  paralleled  with  Sparc  up-regulation  in  bleomycin-induced  tissues.  These  observations 
combining  with  the  results  of  anti-fibrotic  effects  of  Sparc  siRNA  in  fibroblasts  of  the  Tgfbrl 
knock-in  mouse  further  support  a  mutually  regulatory  relationship  between  SPARC  and  TGF-p 
signaling. 

In  summary,  studies  described  here  consistently  demonstrated  that  inhibition  of  Sparc  with  siRNA 
significantly  reduced  collagen  expression  in  both  in  vitro  transgenic  Tgfbrl  fibroblast  model  and  in 
vivo  bleomycin-induced  fibrotic  mouse  model.  This  is  a  first  attempt  to  examine  the  anti-fibrotic 
effects  of  SPARC  inhibition  in  skin  and  lungs  in  vivo.  The  results  obtained  from  these  studies 
provide  favorable  evidence  that  SPARC  may  be  used  as  a  bio-target  for  application  of  anti-fibrosis 
therapies. 
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Figure  legend 

Figure  1.  Comparison  of  gene  expression  between  the  fibroblasts  of  TBR1CA;  Cre-ER  mice 
injected  with  oil  and  4-OHT.  The  expression  level  of  each  gene  in  the  fibroblasts  of  TBR1CA; 
Cre-ER  mice  injected  with  oil  was  normalized  to  1.  Bars  show  the  mean  ±  SD  results  of  analysis  of 
3  independent  experiments  performed  in  triplicate.  *,p  <0.05. 


Figure  2.  Gene  and  protein  expression  in  original  and  siRNA  treated  fibroblasts  from  TBR1CA; 
Cre-ER  mice  injected  with  4-OHT.  (A)  Relative  transcript  levels  of  Col la2,  Ctgf,  and  Sparc.  The 
expression  level  of  each  gene  in  the  fibroblast  lines  with  Non-Targeting  siRNA  (NT  siRNA) 
treatment  was  normalized  to  1.  *,  p  <  0.05.  (B)  Western  blot  analysis  of  type  I  collagen  (Coll),  Ctgf, 
and  Sparc  in  the  fibroblasts  from  constitutive! y  active  Tgfbrl  mice,  with  Ctgf  siRNA  or  Sparc 
siRNA  transfection.  N  =  Non-Targeting  siRNA  treatment;  C  =  Ctgf  siRNA  treatment;  S  =  Sparc 
siRNA  treatment.  (C)  Densitometric  analysis  of  Western  blots  for  protein  level  of  Coll,  Ctgf,  and 
Sparc.  Compared  to  non-targeting  siRNA  treatment,  Ctgf  siRNA  or  Sparc  siRNA  treated  fibroblasts 
showed  significant  reduction  of  Coll  (p  =  0.015  or  0.009  respectively),  Ctgf  (p  =0.002  or  0.056, 
respectively)  and  Sparc  (p  =  0.355  or  0.0004,  respectively). 


Figure  3.  Representative  histological  analysis  of  HE  and  Trichrome  stain  of  mouse  skin  with 
different  treatments  for  4  weeks  in  low  (4  x)  and  high  magnifications  (20  x  ).  A.  Injection  with 
saline  (negative  control)  only;  B.  Injection  with  bleomycin  only;  C:  Injection  with  bleomycin  and 
treatment  with  SPARC  siRNA;  D.  Injection  with  bleomycin  and  treatment  with  Ctgf  siRNA. 


Figure  4.  Representative  histological  features  of  HE  and  Trichrome  stain  of  mouse  lung  samples 
with  different  treatments  intratracheally  in  low  (4  x)  and  high  magnifications  (40  x).  A.  Injection 
with  saline  (negative  control)  only;  B.  Injection  with  bleomycin  only  on  day  0;  C:  Injection  with 
bleomycin  on  day  0  and  SPARC  siRNA  on  days  2,  5,  and  12;  D.  Injection  with  bleomycin  on  day  0 
and  Ctgf  siRNA  on  days  2,  5,  and  12. 


Figure  5.  Gene  expression  in  skin  or  lung  samples  with  different  treatments.  (A)  and  (B)  were  the 
relative  transcript  levels  of  Colla2,  Col3al,  Ctgf,  and  Sparc  in  si  R  N  A -treated  or  untreated 
bleomycin-induced  skins  or  lungs,  respectively.  The  expression  level  of  each  gene  in  the  skin  or 
lung  sample  from  saline  treated  mice  was  normalized  to  1.  *,  p  <  0.05. 


Figure  6.  Collagen  contents  in  skin  (A)  or  lung  (B)  samples  with  different  treatments.  The  collagen 
content  in  the  skin  or  lung  sample  from  saline  treated  mice  was  normalized  to  1.  T reatments:  Sa: 
saline;  B:  bleomycin;  B  +C:  bleomycin  and  Ctgf  siRNAs;  B  +S:  bleomycin  and  Sparc  siRNA.  *,  p 
<0.05. 
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Figure  1. 
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Figure  3. 
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Figure  4. 
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Figure  5. 
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Figure  6. 
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